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Abstract 

The conserved NineTeen protein complex (NTC) is an integral subunit of the spliceosome and required for intron removal 
during pre-mRNA splicing. The complex associates with the spliceosome and participates in the regulation of 
conformational changes of core spliceosomal components, stabilizing RNA-RNA- as well as RNA-protein interactions. In 
addition, the NTC is involved in cell cycle checkpoint control, response to DNA damage, as well as formation and export of 
mRNP-particles. We have identified the Numl protein as the homologue of SPF27, one of NTC core components, in the 
basidiomycetous fungus Ustilago maydis. Numl is required for polarized growth of the fungal hyphae, and, in line with the 
described NTC functions, the numl mutation affects the cell cycle and cell division. The numl deletion influences splicing in 
U. maydis on a global scale, as RNA-Seq analysis revealed increased intron retention rates. Surprisingly, we identified in a 
screen for Numl interacting proteins not only NTC core components as Prp19 and Cefl, but several proteins with putative 
functions during vesicle-mediated transport processes. Among others, Numl interacts with the motor protein Kinl in the 
cytoplasm. Similar phenotypes with respect to filamentous and polar growth, vacuolar morphology, as well as the motility 
of early endosomes corroborate the genetic interaction between Numl and Kinl. Our data implicate a previously 
unidentified connection between a component of the splicing machinery and cytoplasmic transport processes. As the numl 
deletion also affects cytoplasmic mRNA transport, the protein may constitute a novel functional interconnection between 
the two disparate processes of splicing and trafficking. 

Citation: Kellner N, Heimel K, Obhof T, Finkernagel F, Kamper J (2014) The SPF27 Homologue Numl Connects Splicing and Kinesin 1-Dependent Cytoplasmic 
Trafficking in Ustilago maydis. PLoS Genet 10(1): e1004046. doi:10.1371/journal.pgen.1004046 

Editor: Joseph Heitman, Duke University Medical Center, United States of America 

Received May 21, 2013; Accepted October 30, 2013; Published January 2, 2014 

Copyright: © 2014 Kellner et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by grants from the German Research Foundation (DFG) to JK (grants KA 872/6-1 and KA 872/5-1) and the Baden- 
Wurttemberg State Foundation (Landesstiftung Baden-Wurttemberg) to NK. The funders had no role in study design, data collection and analysis, decision to 
publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* E-mail: joerg.kaemper@kit.edu 



Introduction 

The basidiomycete Ustilago maydis is the causal agent of the smut 
disease on corn {^ea mays). During its life cycle, U. maydis displays 
yeast-like, haploid cells that divide by budding, and dikaryotic cells 
that grow as filamentous hyphae. The filamentous stage is initiated 
by fusion of two yeast-like cells (sporidia) and marks the onset of 
the biotrophic stage in which the fungus depends on its host plant 
for propagation. The switch to polarized growth to form the 
elongated filament is indispensable for the successful infection of 
the host plant. 

The regulatory circuits that underlie the dimorphic switch and 
concurrent pathogenic development in U. maydis have been well 
studied within the past years. Similar to other basidiomycetes, a 
heterodimeric complex of two homeodomain transcription 
factors, both encoded by the A-mating type locus, represents 
the major regulatory instance. The two proteins, termed bE and 
bW in U. maydis, form a functional heterodimer, but only if they 
originate from different ^-alleles (e.g. bEl and bW2). Hyphae 
formed upon activation of the i-pathway grow unipolar, but 
only the tip compartment is filled with cytoplasm and contains 
the two genetically distinct nuclei, while the distal part of the 



hypha is composed of empty segments separated by evenly 
spaced retraction septa. Simultaneously with the switch towards 
filamentous growth, cells become arrested in the G2 phase of 
the cell cycle. Only after the penetration of the plant surface, 
this cell cycle block is released, and the "true" filament with 
multiple septated compartments is developed [for recent review, 
see 1]. 

The bE/bW-heterodimer orchestrates a hierarchic, multi- 
layered transcriptional network. Only few /(-responsive genes are 
direct targets of the bE/bW heterodimer. The majority of the 
genes is regulated via the b-dependently induced C2H2 zinc-finger 
transcription factor Rbf 1 . Rbf 1 is the master regulator for several 
other transcription factors that coordinate expression of multiple 
genes associated with cell cycle coordination, morphogenesis and 
pathogenic development. Deletion of rbfl prevents the formation 
of the 6-dependent filaments, and ectopic expression of rbfl is 
sufficient to induce the dimorphic transition. Thus, rbfl is both 
sufficient as well as required for the switch from budding to 
polarized filamentous growth [2]. 

Prerequisite for the growth of the filamentous hyphae is the 
establishment and maintenance of a defined axis of polarity. The 
filaments expand by polar tip growth, which is dependent on long- 
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Author Summary 

In eukaryotic cells, nascent mRNA is processed by splicing 
to remove introns and to join the exon sequences. The 
processed mRNA is then transported out of the nucleus 
and employed by ribosomes to synthesize proteins. 
Splicing is achieved by the spliceosome and associated 
protein complexes, among them the so-called NineTeen 
complex (NTC). We have identified the Numl protein as 
one of the core components of the NTC in the fungus 
Ustilago maydis, and could show that it is required for 
polarized growth of the filamentous fungal cells. Consis- 
tent with the NTC function, cells with a num 7-deletion 
show reduced splicing of mRNA. Moreover, we uncover a 
novel cytoplasmic function of the Numl protein: It 
physically interacts with the microtubule-associated Kine- 
sin 1 motor protein, and phenotypic analyses corroborate 
that both proteins are functionally connected. Our findings 
reveal a yet unidentified role of a global splicing factor 
during intracellular trafficking processes. A possible con- 
nection between these disparate mechanisms presumably 
resides in mRNA-export out of the nucleus and/or the 
transport of mRNA within the cytoplasm. 

distance transport towards the growth cones at the cell apices. This 
directed transport is facilitated by arrays of polarized microtubules 
and a highly conserved set of microtubule-dependent kinesin and 
dynein motor proteins [reviewed in 3]. The cellular cargos that 
rely on microtubule-based transport include endosomes, peroxi- 
somes and nuclei, but also mRNA, which, as shown recently in U. 
maydis, is as well instrumental for establishment and maintenance 
of polarity [4—6] . The latter, however, is not a direct microtubule- 
based cargo, but passively travels on endosomes [7] . In U. maydis, 
transport is mainly mediated through the concerted action of the 
plus-end directed UNC104-like Kinesin- 3 motor protein Kin3, 
which moves endosomes in both directions within the cell along an 
array of antiparallel microtubules, and the minus-end directed 
dynein motor protein Dynl/2 [8—11]. The minus-end directed 
Dyn 1/2 is particularly important at the poles of filaments, where 
unipolar microtubules, with their plus-ends directed to the tip, 
extend from the antiparallel array [9,12]. The conventional 
kinesin motor protein Kinl additionally contributes to the 
establishment of hyphal morphology by transporting Dynein in 
the direction of the microtubule plus-ends within the hyphal apex 
where a loading zone for the retrograde transport processes is 
established [9,12]. Other known cargos for the Kinl motor 
protein include membranous structures; Kinl was previously 
described to be involved in organelle transport [13] and to foster 
the transport of secretory vesicles into the growing tip [14]. More 
recently, transport of the fungal-specific class- 1 7 myosin Mcs 1 was 
shown to dependent directly on Kinl. Mcsl is attached to vesicles 
and contains an N-terminal myosin motor domain fused to a chitin 
synthase region. Anterograde trafficking of Mcs 1 -positive vesicles 
requires both microtubules and filamentous actin and depends on 
Myosin-5 and Kinesin- 1 , which cooperate in delivering vesicles to 
the sites of exocytosis [15-17]. 

In an effort to achieve a better understanding of the complex 
processes required for the establishment and maintenance of the 
dikaryotic hyphae in U. maydis, we employed a candidate approach 
with genes that, based on function or phenotype in other systems, 
were anticipated to be involved in nuclear migration and 
determination of cell polarity. Here, we describe the molecular 
characterization of one of these genes, numl. The numl mutation 
has originally been identified in the basidiomycete Coprinopsis 



cinerea in a screen for mutants affected in the nuclear migration 
during the formation of the heterokaryotic hyphae [18]. 

We show now that Numl is a homologue of SPF27, one of the 
core components of the highly conserved Prpl9/CDC5L or 
NineTeen (NTC) splicing complex [19]. 

Unexpectedly, we identified the conventional Kinesin 1 motor 
protein Kinl [13,20] to interact with Numl. Similar phenotypes 
of Anuml and Akinl hyphae corroborate a functional intercon- 
nection between the two proteins. Our data implicate a previously 
unidentified connection between a component of the splicing 
machinery and cytoplasmic (Kin 1 -dependent) transport processes 
in U. maydis. 

Results 

Numl, a conserved BCAS2/SPF27-like protein, is required 
for polar growth and septum formation 

In the basidiomycete Coprinopsis cinerea, the numl gene was 
identified in a screen for mutants affected in nuclear migration 
during the initial phase of sexual development [18]. Unlike U. 
maydis, C. cinerea grows strictly as a filament, and mating is initiated 
by fusion of two different haploid hyphae harboring compatible 
alleles of the mating type loci. The dikaryotic filament is then 
generated upon migration of "donor"-nuclei into the "acceptor"- 
mycelium. The numl mutation results in strains that are still able to 
donate nuclei, but fail to accept nuclei in compatible mating 
reactions [18]. To examine whether the gene has a conserved 
function in U. maydis with respect to nuclear migration, we set out 
to analyze the potential U. maydis homologue. The U. maydis 
predicted protein Um01682 (MIPS Ustilago maydis Database 
(MUMDB), http:/ /mips. helmholtz-muenchen.de/genre/proj/ 
ustilago/) exhibits 43% identity and 63% similarity with the C. 
cinerea Numl protein, and 33% identity to the human SPF27 
protein (Figure SI). SPF27 homologues are found in all eukaryotic 
clades, with the exception of the Saccharomycetales (Figure S2). 
SPF27 is part of the Prpl9/CDC5L complex, which is an integral 
component of active spliceosomes and required for intron removal 
[21]. Similar to the SPF27 homologues, the U. maydis protein 
harbors a BCAS2 domain (breast carcinoma amplified sequence 2) 
with a so far unknown function and a classical, basically charged 
nuclear localization signal (Figure SI) [22]. 

To determine the function of Numl in U. maydis, the gene was 
deleted in the strains FBI, FB2, SG200 and AB31. FBI (al bl) and 
FB2 (a2 b2) are compatible wild-type strains that, upon mating, 
form filaments on solid media supplemented with charcoal. SG200 
[al mfa'2 bEl b\V2) forms such filaments without a mating partner, 
as the strain harbors a compatible combination of the bEl and 
bW2 genes [23]. For both FB1/FB2 and SG200, cell division upon 
filament formation is stalled in axenic culture, and hyphal 
proliferation depends on the plant host. Polarized growth in 
axenic culture can be best monitored in strain AB31 {al P crgl :bEl I 
P crgl :bW2), which harbors a set of compatible bEl and bW2 genes 
under control of the arabinose-responsive P c , ; y-promoter [24]. 
Upon arabinose-induced expression of bEl/bW2, AB31 switches 
from budding- to polarized growth and forms long filaments 
reminiscent of those formed after fusion of compatible FB 1 and 
FB2 sporidia [2,24]. 

Yeast-like cells of AB3lAnu?nl- and SG2Q0Anuml showed no 
obvious phenotype when grown in axenic culture in complete or 
minimal media (Figure S3). In plant infection assays, the mixture 
of the compatible strains YQXAnuml and FB2Anuml as well as 
SG200Anuml showed reduced disease symptoms compared to the 
respective wild-type strains; both the number of plants with tumors 
as well as the size of tumors decreased (Figure S3). However, 
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microscopic observation of hyphae after plant infection did not 
reveal major morphological differences to wild-type strains (Figure 
S3). The cell-cell fusion of compatible sporidia to generate the 
dikaryotic hyphae, which can be visualized by the appearance of 
"fuzzy" colonies on charcoal containing media, is not altered by 
the numl mutation (Figure S3). However, the phenotype of 
conjugation hyphae upon treatment with synthetic pheromone as 
well as the phenotype of the dikaryotic hyphae showed consider- 
able alterations from wild-type strains. Pheromone treatment of 
Anuml strains led to the formation of branched conjugation tubes 
that are never observed in wild-type strains. Dikaryotic hyphae of 
Anuml strains appeared thicker, branched, and fused sporidia 
displayed bipolar growth; in wild-type strains, however, the fusion 
event is initiated only from a single cell pole (Figures S4 and S5). 
To analyze the phenotype of the hyphae in more detail, we 
employed strain AB3 1 and induced the formation of the filament 
via the expression of bEl/bW2. Similar to the dikaryotic hyphae, 
AB31 Anuml filaments displayed severe alterations. 95% of the 
filaments were significantly shorter than wild-type AB31 filaments 
(average length of 68.9 |J,m compared to 121.7 u,m, t-test: 
p = 6.35x10 4b ). Filament formation upon ^-induction in AB31 
is initiated at one of the cell poles of the sporidium; the resulting 
hyphae grow unbranched and strictly unipolar (Figure 1A). In 
ABU Anuml, hyphae generally exhibited an irregular and curved 
morphology (Figure IB, Table 1). In 18% of the sporidia, hyphae 
initiated at both cell poles (bipolar, Figure 1C, Table 1), and 4.4% 
exhibited a branched morphology (N = 461) (Figure ID, Table 1). 
In addition, 21.2% of the AB3lAnu?nl filaments had delocalized 
septa, and 2.4% of the hyphae generated empty compartments at 
the cell tips (Figure IE and F, Table 1). Septa in AB3lAnuml could 
be visualized by different chitin-stains, as Calcofluor [25], Congo 
Red [26] or fluorescein-conjugated wheat germ agglutinin (WGA) 
[27] (Figures IE and S6; see also Figure S6 legend). However, in 
contrast to septa in the dikaryotic hyphae, septa in ABSlAnuml 
were not always visible in the DIC channel, and staining with 
fluorescein-WGA revealed extensive chitin accumulations or ring- 
like chitin structures throughout the hyphae, indicating that these 
septa or septa-like structures may be structurally different from 
wild-type septa. 

To ensure that the observed phenotype is due to the numl 
deletion, a construct harboring a numl:eGFP fusion gene under 
control of the arabinose inducible P^y-promoter was introduced 
in single copy into the cbx locus. Strain UNK220 (AB31 Anuml 
P crg i:numl:egfp) was grown in arabinose-containing media to induce 
both the bEl/bW2 heterodimer as well as the numl :eGFP gene; the 
induced filaments were indistinguishable from that of the wild-type 
AB31 strain (Figure S7). 

Numl is a component of the highly conserved 
spliceosome-associated Prp19/CDC5L-complex 

Numl has striking similarities to the human BCAS2/SPF27 
protein that has been identified as one of the four core components 
(Prpl9, PLR1, CDC5L and SPF27) of the Prpl9/CDC5L 
spliceosome-associated complex, commonly termed NineTeen 
Complex (NTC) (Figure 2A). To assess whether Numl is a 
structural component of the putative NTC complex in U. maydis, 
we employed a directed yeast two-hybrid approach using the full- 
length Numl protein as bait as well as full-length and C- and N- 
terminal versions of the Prpl9 (Uml0027) and CDC5L-homo- 
logues (termed Cefl, Um04411) as prey, respectively. Consistent 
with the previously identified interactions in S. pombe [28], Numl 
interacts with the Prpl9 N-terminus as well as the Cefl C- 
terminus (Figure 2B). The interaction of Numl with Prpl9 and 
Cefl was verified by in vivo co-immunoprecipitation. numl was 
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Figure 1. Numl is required for polar growth and septum 
formation. Filament formation of AB31 wild-type and AB31 Anuml 
deletion strains was monitored 1 2-1 4 hours upon induction of the bE1 / 
bW2 heterodimer. (A) AB31 wild-type cells grow as long, straight 
filaments; only the tip compartment is filled with cytoplasm, separated 
from highly vacuolated distal compartments by a septum (arrowhead). 
(B) AB31zlnum7 hyphae grow shorter and more curved. In AB31zfnum!, 
bipolar (C) and branched (D) hyphae are observed frequently. (E) Septa 
of AB31zlnum7 are often formed within the compartment filled with 
cytoplasm (arrowheads denote delocalized septum; in lower panel, 
Calcofluor White staining was used to visualize the septum). (F) In 
approx. 2.5% of hyphae delocalized septa (arrowhead) lead to empty tip 
compartments. Scale bars: 10 urn. 
doi:10.1371/journal.pgen.1004046.g001 

replaced by an eGFP-tagged version of the gene in strain AB3 1 
(a2, P crg :bW2, bEl), and additionally prpl9 or cefl were exchanged 
with HA-tagged versions, resulting in strains UNK200 (AB31 
numl:3xeGFP:hyg R , prpl9:3xHA:cb}f) and UM08 (AB31 nu- 
ml :3 xeGFP:hyg R , cefl :3xHA:cbx R ), respectively. All genes were 
introduced via homologous recombination into their native loci to 



Table 1. Quantification of phenotypes in AB31 and 
AB31 Anuml 12 hours after the induction of the bE/bW- 
heterodimer. 







AB31* 


AB31^//7um7* 


long and straight hyphae 


92% ±0.0 10 


4.7%±0.029 


irregular growing hyphae 
(short, curvy) 


8%±0.015 


49.13%±0.029 


bipolar 


0% 


18.13%±0.026 


branched 


0% 


4.4%±0.01 1 


delocalized septa 


0% 


21.2%±0.005 


empty compartments at hyphal 
tips 


0% 


2.43%±0.039 




N = 375 


N = 461 


0 length of hyphae** 


121.7 um±15.63 |im 


68.9 |im± 15.37 urn 




N = 75 


N = 75 



*Mean values from three independent biological replicates are shown. 

**t-test p = 6.35x10~ 46 . 

doi:1 0.1 371 /journal.pgen.1 004046.t001 
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Figure 2. Numl interacts with components of the Prp19/Cdc5L complex. (A) Schematic representation of the protein-protein interactions 
(blue arrows) within the human hPrp19/CDC5L-complex [21]. Proteins depicted in blue represent the core complex, additional components (yellow) 
are associated with individual proteins. Green arrows indicate U. maydis protein interactions identified in this study. (B) The Prp19 N-terminus and the 
Cefl C-terminus interact with the full length Num1 protein in a yeast two-hybrid analysis. p53/SV40 T-antigen and Lamin/SV40 T-antigen (Clontech) 
serve as positive and negative control, respectively. (C) In vivo co-immunoprecipitation of HA-tagged Prp1 9 and Cefl and 3eGFP-tagged Num1 . rx-HA 
coupled agarose beads were used to precipitate Prp19 and Cefl, respectively. The Numl protein was detected on a Western blot using an a-GFP 
antibody. The negative control shows no unspecific binding of the Numl protein to ot-HA agarose. (D) Subcellular localization of the Num1:3eGFP 
fusion protein. In addition to the nuclear localization of Numl (enclosed by dashed circle), a cytoplasmic distribution to distinct foci can be observed 
(arrowhead). Scale bar: 5 am. (E, F) Num1:3eGFP as well as Cef1:RFP and Prp19:RFP were co-expressed in AB31 under control of their native 
promoters. Numl co-localizes with both Cefl and Prp19 in the nuclei. Additional cytoplasmic Num1:3eGFP foci are indicated by arrowheads. Scale 
bars: 10 urn. 

doi:10.1371/journal.pgen.1004046.g002 



ensure native protein levels. All fusion proteins were shown to be 
functional (Supporting Text SI, Figure S8). By use of an anti-HA- 
antibody, Numl:3eGFP was specifically co-precipitated with both 
Prpl9:HA and Cefl:HA in protein extracts from strains UNK200 
and UM08 (Figure 2C). The interaction of Numl with Prpl9 and 
Cefl was further corroborated by their localization. Numl:3eGFP 
and Prpl9:RFP or Cefl:RFP were co-expressed in strain AB31, 
and all proteins co-localized in the nuclei (Figure 2E and F). While 
for Prpl9:RFP and CefhRFP the fluorescence signal was exclu- 
sively localized in the nucleus, Numl:3eGFP was also detected in 
the cytoplasm (Figure 2D). Transformation of AB3\ Anuml with a 
numheGFP derivative harboring a mutated nuclear localization 
sequence (KKRK to AAAA) did not lead to complementation of the 
Anuml phenotype, and the fusion protein localized predominantly 
in the cytoplasm (Figure S7), which implies that nuclear localization 
of Numl is required for function of the protein. 

In addition to their function in mRNA splicing, NTC and 
NTC-associated proteins are known to be involved in several other 
cellular processes. In S. pombe and S. cerevisiae, Cdc5 and Ceflp, 
respectively, were originally identified as cell division cycle 
mutants that are required for G2 progression and mitotic entry 
[29-31], and Prpl9p has been isolated in S. cerevisiae in a screen for 
mutants conferring sensitivity to DNA damaging agents [32]. It 
has been proposed that Prpl9p homologues play a direct role in 
DNA damage repair that is conserved from yeast to human [33— 
38]. 



To define whether Numl acts in line with previously known 
functions of the NTC-complex, we investigated cell cycle 
regulation and the response to DNA damage in AB31 Anuml. In 
AB31, the filaments generated after induction of the bEl/bW2 
heterodimer are arrested in the G2 phase of the cell cycle, 
harboring a single haploid nucleus that is positioned in the tip 
compartment [39]. In filaments of ABZlAnuml, however, devia- 
tions from this pattern became evident. 18% of the filaments 
(N= 100) exhibited two (or more) nuclei within a single hyphal 
compartment; additionally, compartments devoid of any nuclei were 
generated, most likely due to the aberrant formation of septa 
(Figure 3A). In accordance with this observation, FACS analysis of 
AB3 1 cells revealed that upon prolonged induction of the bEl /bW2- 
heterodimer cells with a 2C DNA content accumulated, whereas in 
AB31 Anmnl, an increased number of cells with both 1C and 2C 
DNA content was observed (Figure S9), indicative for a deregulated 
cell cycle. Taken together, our results implicate that the numl 
deletion affects both cell cycle regulation as well as cell division. 

Consistent with the function of the NTC in response to DNA 
damage, we found that Anuml cells exhibited increased sensitivity 
to UV irradiation and DNA-damaging agents. The survival rate 
upon UV-treatment was reduced approximately by a factor of 1 0 
in AB3\ Anuml (Figure 3B). Similarly, treatment with phleomycin, 
causing double-strand breaks in DNA [40] or hydroxyurea, which 
triggers stalling of the replication fork during DNA synthesis [41], 
led to reduced survival of Anuml cells (Figure 3B). 
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Figure 3. Numl influences the cell cycle and cell division as 
well as the cellular response to DNA damage. (A) Induction of the 
bE1/bW2 heterodimer in wild-type cells (AB31) induces a G2 cell cycle 
arrest, resulting in hyphae with a single nucleus positioned in the tip 
compartment. In AB31 Anuml , 18% of the hyphae contain more than 
one nucleus, usually separated by a septum or a septum-like chitin 
structure, which is not always visible in the DIC-channel (see also Figure 
S6). In some cases, delocalized septa give rise to compartments that are 
devoid of any nuclei (bottom panel). Filaments were stained with DAPI 
and Congo Red, which effectively stains chitin in fungal cell walls [26], 
to visualize nuclei and septa within the same hypha. Scale bars: 10 am. 
(B) Sensitivity of an FB2Anum1 deletion strain to DNA-damage in 
comparison to wild-type cells. Serial 10-fold dilutions were spotted on 
complete medium containing the indicated DNA-damage inducing 
supplements, or subsequently irradiated with 300 J/m 2 UV light 
(254 nm), respectively. Pictures were taken 2-3 days after incubation 
at 28°C. 

doi:10.1371/journal.pgen.1004046.g003 



Num1 affects splicing of mRNA on a global scale 

To assess whether the Anuml mutation affects splicing, we 
employed next-generation RNA sequencing (RNA-Seq) to deter- 
mine the transcriptome of strains AB31 and AB3lAnuml. bEl I 
bW2 expression was induced for eight hours; at this time point, in 
both strains a comparable quantity of filamentous cells was 
observed. The expression levels of b W s of the central regulator rbfl 
(a direct A-target gene) and oinuml were monitored by quantitative 
RT-PCR analysis to ensure similar expression levels (Figure S10). 
For both AB31 and ABZlAnuml, RNA-Seq experiments from 
three independent biological samples were performed, and 
sequences were compared to the manually curated U. majdis gene 
models of the MIPS Ustilago maydis Database (MUMDB). Splicing 
efficiency was calculated from the normalized read number within 
introns in relation to the normalized read number of exon regions 
of a gene (for details, see Materials and Methods). Based on 2142 
intron sequences, introns of AB31 genes were spliced with an 
average efficiency of 0.93 (±0.13); splicing efficiency in AB31z(- 
numl was significandy lower (factor — 1.44, t-test p = 6.22*10 -111 ), 
with an average of 0.80 (±0.23). The splicing defect observed by 
RNA-Seq was verified for a small set of genes by qRT-PCR 
(Figure Sll). In general, we noticed a wide variation in splicing 
efficiency of different introns in response to the 7zwm/-deletion. For 
80.2% of the introns (1717 of 2142), the intron retention rate 
(fraction of RNA in intronic sequences in relation to exonic 
sequences) in AB'il Anuml was at least twice as high as in the wild- 
type strain, indicating that Numl is required for splicing in 
general. We noticed, however, that for few introns (83, allowing a 
10% threshold) splicing was not affected, while other introns were 
barely spliced at all in AWil Anuml. Differences in intron retention 
rates were specific for discrete introns, and not for particular 
transcripts, as we observed that in genes with more than one intron 
separate introns were spliced to different extent (Table SI, 
examples are shown in Figure SI 2). 

It has been proposed that the function of the NTC with regards 
to cell cycle regulation or DNA damage repair may result from the 
decreased splicing efficiency of intron-containing genes involved in 
recombination, repair, cell cycle or chromosome segregation [42] . 
We did not observe an enrichment of genes with high intron 
retention rates in any "category for functional annotation" 
(FunCat). However, inspection of genes in the FunCats "cell 
cycle" and "DNA repair" (Table S2) revealed that for individual 
genes involved in DNA repair splicing is impaired, as, for example, 
for rad51 (um()3290), for the gene for the repair exonuclease Reel 
(uml 1016), or for mms2 (uml0097), involved in the error-free post- 
replication repair. Similarly, splicing of genes involved in cell cycle 
control, as cksl (um03210), encoding a regulatory subunit of cyclin- 
dependent kinases, or the gene for the cyclin dependent kinase 1 
[cdkl, urn 10705) is affected in AB3lAnuml (Table S2, Figure S13). 

We have shown before that cell cycle and hyphal development 
are controlled by a cascade of transcription factors [2]. Intrigu- 
ingly, although most U. maydis genes contain no introns, both bW 
and bE cany single introns; rbfl harbors even four introns. It could 
be well possible that the observed phenotype of the Anuml strains 
is caused by inefficiently spliced genes of the master regulators for 
filamentous growth and pathogenic development. In AB31, the 
b W2- allele is expressed without the intron sequence [24], but we 
observed that bEl is indeed spliced less efficiently in AB.3lAnHmi 
(Table S3). However, all direct b-target genes were induced to a 
similar extent in AB31 Anuml compared to AB31 (Table S4), 
which implies that expression of the bE and bW genes in the Anuml 
background is sufficient to provide full functionality of the bE/bW 
heterodimer. 
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In accordance with the unaltered functionality of bE/bW, the 
expression of the direct target gene rbfl, was not altered in 
AB'ilAnuml (Figure 4C); however, all four introns were spliced less 
efficiendy . Under the assumption that splicing of the introns occurs 
independendy from each other, only in 45 % of the mature mRNA 
in ABil Anuml all introns were spliced correctly, compared to 
80% in AB31 wild-type cells (Figure 4A, Table S3). The splicing 
defect in rbfl was validated via RT-PCR and qRT-PCR analyses 
(Figure SI 4), showing additional unspliced variants and increased 
intron retention, respectively, in AB'ilAnuml. Consistently, we 
observed that the abundance of the Rbfl protein in ABHAnuml is 
reduced to about 30% of wild-type level (Figure 4B), although 
qRT-PCR analysis showed that the total level of the rbfl -transcript 
was not altered (Figure 4C). 

Rbfl controls the expression of bizl, hdpl and hdp2, encoding 
transcription factors that are collectively involved in cell cycle 
regulation, filamentous growth and pathogenic development [43, 
Scherer, Vranes, Pothiratana and Kamper, unpublished]. In 
accordance with the reduced Rbfl levels, RNA-Seq analysis 
revealed that bizl as well as hdp2 expression was significandy 
down-regulated in AB'ilAnuml , which was verified by qRT-PCR 
(Figures 4D and SI 5, Table S4). hdpl splicing efficiency was 
reduced in AB31 Anuml from 0.83 to 0.14 (Table S3); interestingly, 
hdpl RNA-levels were increased about 2-fold, which may argue for 
a self-regulation of the gene (Figures 4D and SI 5, Table S5). In 
total, we observed an at least two-fold down-regulation for 228 
genes in AB'MAnuml; intriguingly, 89 of these genes overlap with a 
set of 351 genes that has been identified as down-regulated in 
response to a deletion of rbfl in AB31 (Figure 4E, Table S5, [2]). 

As rbfl is one of the major regulators implicated in filamentous 
growth and cell cycle regulation, it was tempting to speculate that 
the reduced splicing efficiency of the gene may be the cause for the 
observed phenotype in zinami-mutants. To test this hypothesis, we 
exchanged the rbfl gene with an intron-free derivative in strains 
AB3 1 and AB'ilAnuml. However, we observed that the intron-free 
copy of rbfl is not capable to substitute for the native gene; despite 
wild-type-like expression levels the strains did not grow filamen- 
tous, nor were the known rbfl -target genes induced (Figure SI 6). 
Obviously, introns are essential for the expression of the Rbfl- 
protein. As a consequence, expression of the rbfl cDNA did not 
rescue the Anuml phenotype either. 

In summary, we conclude that Num 1 is a structural component 
of the conserved Prpl9/CDC5L complex and fulfills common 
functions that have previously been associated with this complex, 
as cell cycle control, cellular response to DNA damage, and 
splicing. Our data indicate that the observed phenotype in U. 
maydis hyphae may be at least in part the result of the splicing 
defect of genes involved in cell cycle control, DNA damage repair 
and of central components of the regulatory cascade controlling 
filamentous growth and pathogenic development. 

Num1 interacts with the Kinl motor protein 

To further elaborate the molecular function of the Numl 
protein, we employed the yeast two-hybrid system. Proteins 
identified as potential interactors for Numl are shown in Table 
S6. Consistent with a putative function of Numl within the 
Prpl9/CDC5L complex, the CDC5L homologue Cefl was 
isolated as an interacting protein. 

Unexpectedly, we identified various interactions with proteins 
involved in cellular transport processes. Uml0158 is an adaptin- 
like protein from clathrin/coatomer-adapter complexes 
(IPR002553) mediating endocytic protein transport between ER 
and Golgi. Um03539 and Umll510 contain a BAR-domain 
(IPR004148) and a BROl -domain (IPR004328), respectively, and 



are thus likely involved in intracellular vesicular transport 
processes and protein targeting to endosomes or the vacuole. 
The conventional kinesin motor protein Kinl [formerly named 
Kin2, 20] functions in long distance trafficking within fungal 
hyphae as it transports Dynein along the microtubule cytoskeleton 
towards their plus-ends directed to the hyphal tips [9, 1 4] . We were 
not able to generate deletion strains for uml0158, indicating that 
the respective protein might be of essential function. Deletion 
analysis for um03539 and umll510 in strain AB31 revealed no 
obvious phenotype with respect to polar growth and vesicular 
movement (data not shown), indicating that the respective proteins 
have no function with respect to Num 1 -associated processes. For 
kinl, however, it had already been demonstrated that the gene has 
impact on hyphal morphology [14, see below]. 

To confirm the interaction of Num 1 with Kin 1 , co-immuno- 
precipitation analyses were carried out. A full-length Num 1 :Myc- 
tagged protein and a Kinl:HA-tagged protein fragment encom- 
passing amino acids 650-968 were generated by coupled in vitro 
transcription/translation, and Numl:Myc was specifically co- 
precipitated from mixtures of Kinl:HA with an anti-HA antibody 
(Figure 5 A). By use of the yeast two-hybrid system we mapped the 
interaction domain of Numl (encompassing amino acids 67 to 
147) with the C-terminal Kinl domain used in the co-immuno- 
precipitation, which provides additional evidence for a specific 
interaction with the Kinl motor protein (Figure SI 7). The 
interaction was further verified in vivo, using U. maydis strain 
UNK197 that harbors numl:3eGFP and kinhHA fusion genes 
integrated into their native loci by homologous recombination, 
respectively. By use of dithiobis [succinimidyl] -propionate (DSP) to 
crosslink proteins in extracts prior to immunoprecipitation, the 
Numl:3eGFP protein was efficiently co-precipitated with the 
Kinl: HA protein using anti-HA antibodies (Figure 5B). Without 
cross-linking of the proteins, only a faint signal was detectable (not 
shown), suggesting that the interaction is only transient. Taken 
together, we conclude from our data that Numl interacts with the 
motor protein Kinl. 

Numl and the motor protein Kinl are functionally 
connected 

In addition to the physical interaction, several lines of evidence 
argue for a functional connection between Numl and Kinl. 
Induced expression of bEl/bW2 in a Akinl strain leads to short 
curved hyphae with delocalized septa and the formation of bipolar 
filaments [14]. Similarly, the induction of hyphal growth in 
AB'ilAnuml leads to irregular or bipolar filaments with an 
aberrant insertion of septa (Figure 1, Table 1). In strains deleted 
for kinl the number of vacuoles is increased, while size is reduced 
[13]. We observed a similar phenotype in strains SG200Anuml 
(Figure 6A). 

To analyze the functional relation between the two proteins 
further, we generated strains deleted for both genes. Filaments 
induced in ABil Anuml / Akinl displayed a more severe phenotype 
with respect to polarized growth, as hyphae exhibited a very 
irregular, swollen morphology, and also the number of branched 
hyphae and hyphae with an altered septation pattern was 
significandy higher than in AB'ilAnuml or ABilAkinl (Figure 
S18). Infection experiments with SG200zfnam7, SG2Q()Akinl and 
SG200 A numl / Akinl deletion strains revealed that the virulence of 
the double deletion strain is significantly reduced when compared 
to the respective single deletion strains (Figure SI 8). The more 
severe phenotype excludes that both genes are in pure epistatic 
interaction; besides their function within a common complex, both 
proteins must fulfill functions that are independent from each 
other. 
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Figure 4. Reduced splicing efficiency of the rbfl-gene leads to impaired function of the Rbfl master regulator. (A) Plot depicting 
splicing efficiency of the rbfl-gene. Plotted are the FPKM values (fragments per kilobase of sequence per million fragments mapped) across the 
genomic region indicated (coordinates in nucleotides) of three independent RNA-Seq experiments for AB31 wild-type (blue lines) and AB31 Anuml 
(red lines), respectively. Exons (E) and introns (I) are indicated. All four introns show increased intron retention rates in AB31 Anuml. (B) Western 
analysis showing abundance of Rbf1:3xHA and a-tubulin (loading control) from AB31 wild-type and Anum7-deletion strains. In AB31z/num7, Rbfl is 
reduced to 30% of wild-type-level (Quantification: ImageJ [109]). (C, D) Gene expression analyses of b- and rbfl -genes (C) as well as /-far"; -target genes 
(D) using qRT-PCR. RNA samples were isolated from strains AB31 and AB31 Anuml eight hours after induction of the bE1/bW2-heterodimer. Gene 
expression is shown relative to the highest expression value, using actin and elF2b for normalization. Shown are the mean values of three biological 
and two technical replicates. Error bars represent the SD. (E) Venn diagram depicting the total number of genes repressed in AB31 Anuml and 
AB3'\Arbf1. * RNA-Seq analysis, this study; ** Microarray analysis conducted five hours after induction of the bE1/bW2-heterodimer [2]. 
doi:10.1371/journal.pgen.1004046.g004 



Num1 and Kinl are both required for trafficking of early 
endosomes 

The deletion of kinl results in a loss of rapid bi-directional 
trafficking of early endosomes. Early endosomes visualized by a 
GFP-fusion with the endosomal t-SNARE Yupl (YuplxGFP) [4] 
accumulate at the hyphal tips of Ami-deletion mutants [9], while in 
the wild-type situation YuplxGFP labeled vesicles shuttle bi- 
directionally throughout the hyphae along microtubule tracks [4]. 
The anterograde transport of early endosomes is mediated via the 
motor protein Kin3 (directed to the microtubule plus-ends at the 
hyphal tip), whereas the retrograde transport is accomplished by 



Dyneinl/2 [8,10,11]. Since the Kinl motor protein transports 
Dynein toward the hyphal tips, the accumulation of early endosomes 
at the tip of Ndnl hyphae can be explained indirectly by the failure of 
the Kin 1 -dependent anterograde transport of Dynein [9] . 

The deletion of numl leads to a comparable scenario in the 
distribution of early endosomes, which appeared to cluster at 
hyphal tips and the distal pole of the hyphae as well as around the 
delocalized septa Frequently found in A man 1 -hyphae (Figure 6B). 
Endosomal motility was monitored in strain AB33, in which the 
&-genes are expressed under control of the nitrogen-responsive 
-Pjiarrpromoter [44]. Filamentous growth in AB33 is repressed by 
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Figure 5. Numl physically interacts with the Kinl motor 
protein. (A) In vitro expression and co-immunoprecipitation of Myc- 
tagged Num1 and HA-tagged Kin1. Proteins were labeled with biotin 
during in vitro synthesis, detected with streptavidin-conjugated alkaline 
phosphatase, and visualized with NBT/BCIP. In the negative control 
(Num1:Myc), no unspecific binding of Myc-tagged Num1 to the HA- 
antibody was detected. In contrast, after co-incubation with HA-tagged 
Kinl, the Myc-tagged Num1 protein was co-immunoprecipitated 
(indicated by asterisk), demonstrating that the Numl protein physically 
interacts with Kinl in vitro. (B) In vivo co-immunoprecipitation of U. 
maydis strains expressing HA-tagged Kinl and 3eGFP-tagged Numl 
under control of their endogenous promoters. After crosslinking of 
proteins using DSP, a-HA coupled agarose beads were used to 
precipitate Kinl . The Numl protein was detected on Western blots 
using an a-GFP antibody. As negative control, a strain expressing only 
Kin1:3xHA was used. 
doi:1 0.1 371 /journal.pgen.1 004046.g005 

ammonium and can be elicited in nitrate-containing media. In 
AB33Anuml filaments, long-distance movement of early endo- 
somes was constrained, and velocity of the vesicles was reduced to 
a speed of 1.14 |im/sec in AB33Anuml, in contrast to 1.46 um/sec 
in wild-type hyphae (t-test p = 0.005). In AB33Anuml filaments 
with aberrant septa, long-distance trafficking of endosomes was 
entirely abolished, and only residual movement of vesicles was 
observed (Table 2, Videos SI and S2). Next we addressed the 
question whether the altered movement and positioning of early 
endosomes in Anuml strains may be caused by a distortion of 
cytoskeletal elements. Neither the arrangement or abundance of 
microtubules (visualized by an a-tubulin mCherry fusion, 
Tub 1 :mCherry) nor of actin patches (visualized by LifeAct: YFP 
[45]) was altered in the Anuml mutant strain (Figure S19). To 
analyze the orientation of the microtubule cytoskeleton, we used 
strains expressing an a-tubulin GFP fusion [14] together with the 
microtubule plus-end marker Peb 1 , fused to RFP, which localizes 
to growing microtubule plus-ends [46]. As previously described 
[46], in wild-type interphase sporidia the majority of microtubule 
plus-ends is oriented towards the budding cell and the opposing 
cell pole of the mother cell; in Anuml sporidia, the orientation of 
the microtubule cytoskeleton was not considerably altered (Figure 
S20 A,B). Similarly, in hyphae the majority of microtubule plus- 
ends is oriented towards the cell poles, i.e. hyphal tips and basal 
retraction septa [9]. In A numl -mutants, we observed microtubule 
plus-ends also in close vicinity to the delocalized septa (Figure S20 
C-E). Apparendy, endosomes accumulate at microtubule plus- 
ends in mm 1 -deletion mutants, supporting the Kin 1 -dependent 
effect on endosomal transport. In some hyphae of the numl- 
deletion strain we observed the appearance of mitotic spindles, 
which supports the defect in cell cycle regulation. 



To determine whether the effect on endosomal trafficking in 
AB33Anuml strains is attributed to an aberrant anterograde 
transport of Dynein towards the hyphal apex, as it was suggested 
for Akinl mutants, we visualized the subcellular distribution of 
Dynein using an N-terminal 3eGFP:Dyn2 fusion protein [9]. In 
wild-type filaments, 3eGFP:Dyn2 fusion proteins form comet-like 
motile structures that accumulate at microtubule plus-ends in 
hyphal tips, where they reverse direction and move into retrograde 
direction [9]. In AB33Anuml filaments, however, the number of 
cells harboring an apical accumulation of 3eGFP:Dyn2 was 
significandy reduced (t-test p = 0.002); concomitantly, in most cells 
the 3eGFP:Dyn2 signal in the cytoplasm was higher than in wild- 
type filaments (Figure 6C). 

As Kin 1 affects motility of early endosomes only indirectly (via 
the anterograde transport of Dynein), we set out to analyze the 
movement of a direct Kinl -cargo. The myosin chitin synthase 
Mcsl [15] localizes to Kinl -dependently transported vesicles [17]. 
It has been described before that 3eGFP:Mcsl fusion proteins 
localize to the poles of growing sporidia, where it is secreted and 
participates in the synthesis of chitin [15]. A num 1 -mutants showed 
reduced Mcsl -accumulations at the growth region in budding 
cells, which coincides with the distribution of Mcsl in fan /-deletion 
mutants [17] (Figure S21 A,B). In AB33 filaments, 3eGFP:Mcsl 
localizes to distinct foci close to the cell membrane and forms a 
gradient towards the growth zone within the hyphal apex or 
localizes to the basal retraction septum, i.e. to sites of chitin 
synthesis. In contrast, in AB33Anuml hyphae no tip-directed 
gradient was obvious and less foci were observed at the cell 
membrane (Figure S21 C,D), indicating an impaired delivery of 
Mcs 1 -containing vesicles. In many cases, Mcsl was not observed 
at delocalized septa in Anuml-mutant hyphae. In summary, U. 
maydis strains deleted for either numl or kinl show similar 
phenotypes with respect to hyphal morphology and septation, 
formation of vacuoles, endosomal trafficking as well as the 
localization of the Kin 1 -cargos Dynein and Mcs 1 , which strongly 
corroborates a functional interrelation of Num 1 and Kin 1 . 

Numl is required for transport of cytoplasmic mRNA 
particles 

Recently, it has been shown that the RNA-binding protein 
Rrm4 is associated with early endosomes [7] . Rrm4 is required for 
cytoplasmic transport of mRNAs in U. maydis filaments [47]. We 
deleted numl in an AB33 strain expressing an Rrm4:eGFP fusion 
protein. In the wild-type background vesicles associated with 
eGFP-labeled Rrm4 shuttled within the hyphae to the distal and 
proximal ends. In contrast, movement in Anuml was reduced, and 
similar to Yupl -labeled early endosomes, Rrm4-containing 
particles frequendy accumulated in clusters at hyphal tips or 
around delocalized septa (Figure S22, Videos S3, S4), indicating 
that numl is required for cytoplasmic mRNA transport. 

Numl has a functional homologue in the filamentous 
ascomycete Aspergillus nidulans 

To address the function of Numl in filamentous fungi in 
general, we deleted a gene encoding a protein with high 
similarities to Numl (AN4244, 49% similarity, 30% identity) in 
the ascomycete Aspergillus nidulans (strain TN02A3). The deletion 
mutants displayed smaller colonies with a reduced rate of 
conidiospore formation. Germination of spores was drastically 
reduced in the AAM4244 mutant strains: after eight hours of 
incubation at 37°C, only 5% of the spores germinated, compared 
to 98% in the wild-type strain. Similar to the phenotype observed 
in U. maydis, hyphae from germinating AAM4244 conidia were 
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Figure 6. Numl and Kin1 are functionally connected. (A) The num1 deletion leads to aberrant vacuole morphology. Vacuoles were visualized 
by CellTracker Blue (7-amino-4-chloromethyl-coumarin, CMAC) staining. Both numl and fan 7 deletion strains contain more, but smaller vacuoles in 
comparison to wild-type sporidia. Scale bars: 5 jim. Right panel: Quantification of the average number of vacuoles of strains indicated. N represents 
the number of individual cells analyzed. Mean values of three independent experiments are shown. (B) The numl deletion leads to an aberrant 
distribution of early endosomes. For visualization of early endosomes, a Yup1:eGFP fusion protein [9] was expressed under control of the 
constitutively active P oter promoter in strains AB33 and AB33Anum1, which harbor b-genes under control of the nitrate-inducible P nar ,-promoter. 
Endosomes were analyzed 14 hours after bE J/bW2-induction. Septa (arrowhead) were visualized by Calcofluor White staining. Yupl -labelled early 
endosomes accumulate at basal and apical parts of the hyphae as well as around delocalized septa in AB33zfnum7 in contrast to wild-type hyphae. 
Scale bars: 10 p.m. Movies detecting Yup1 :eGFP-fluorescence of the depicted hyphae are shown in Videos SI and S2. Shown below are the 
corresponding kymographs (time and distance as indicated) reflecting the motility of early endosomes. In AB33 hyphae, early endosomes travel 
continuously towards the hyphal tip (anterograde) and reverse direction to travel back towards the septum (retrograde). Pausing organelles are 
reflected by a vertical line. In contrast, endosomal motility is drastically reduced in AB33z)num7 hyphae containing additional septa, where the 
majority of the Yup1:eGFP-signal is non-motile. An arrowhead indicates the septum. (C) The localization of Dynein is dependent on Numl. In AB33 
filaments 3eGFP:Dyn2 fusion proteins form comet-like motile structures that accumulate at the hyphal tip and from there move into retrograde 
direction [9], In filaments of AB33zlnum/ strains the apical accumulation is significantly reduced (p = 0.002); instead a stronger signal is observed in 
the cytoplasm. Dynein localization was monitored 14 hours after induction of the bE1/bW2 heterodimer. Scale bars: 10 |im. Right panel: 
Quantification of the apical accumulation of 3eGFP:Dyn2 fusion proteins. N represents the number of hyphae analyzed. Mean values of three 
independent experiments are shown. 
doi:10.1371/journal.pgen.1004046.g006 
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Table 2. Reduced endosomal motility in num 7-deletion strains. 
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shorter and thicker than those of the wild-type strain. Staining with 
CellTracker Blue revealed that in the mutant the number of 
vacuoles was increased, while size was reduced (Figure 7A). In 
contrast to V. maydis, A. nidulans hyphae have cell compartments 
with multiple nuclei. In wild-type hyphae, these nuclei are evenly 
distributed. In AAN4244, however, nuclei were often found in 
clusters (Figure 7B). A similar phenotype has been observed earlier 
for strains deleted for KinA, the homologue of the U. maydis Kin 1 
[48] . To test a potential functional relation between AN4244 and 
KinA in A. nidulans, we analyzed, in analogy to U. maydis, the 
distribution of Dynein (NudA in A. nidulans). We deleted AN4244 in 
strain LZ12, which harbors a GFP:NudA fusion protein. Similar 
as described for the U. maydis hyphae, GFP:NudA was predom- 
inantly localized at the tip of hyphae in strain LZ12, while the 
apical accumulation was reduced in LZ 12AAN4244 (Figure 7C). 
Taken together, our data suggest that the A. nidulans AN4244 has a 
conserved function as the U. maydis protein with respect to 
polarized growth and intracellular transport processes. 

Discussion 

In this study we have characterized the function of the Num 1 
protein in U. maydis. We show that Numl is a functional 
homologue to the human SPF27 protein, which is an integral 
component of the conserved spliceosome-associated Prpl9/ 
CDC5L-complex (NTC). Deletion of the numl gene affects polar 
and filamentous hyphal growth, cell cycle control as well as DNA 
repair. Consistent with its function as component of the NTC, 
deletion of numl results in a genome-wide reduced splicing- 
activity, which could explain the observed pleiotropic phenotype. 
In addition we demonstrate a previously unknown function of an 
NTC-component in the cytoplasm, where Numl interacts with the 
Kinesin 1 motor protein. This functional interconnection of Numl 
with elements of the cytoskeleton implies additional, yet unknown 
functions during intracellular trafficking processes. 

Numl is a conserved component of the spliceosome- 
associated Prp19/CDC5L-complex (NTC) 

Numl shows significant similarities to SPF27, and its interaction 
with Prpl9 and Cefl is in accordance with the interactions within 
the NTC complex demonstrated in S. pombe [28], humans [21] and 
plants [49]. 

The NTC was initially identified in S. cerevisiae as a complex 
containing Prpl9 and at least seven other components [50], most 
of which are conserved from yeast to human. The NTC is a non- 
snRNP (small nuclear ribonucleoprotein particle) multi-protein 
complex, and as a subunit of the spliceosome it is required for 
intron removal during pre-mRNA splicing. The complex associ- 
ates with the spliceosome during or after dissociation of U4 snRNP 
and participates in the regulation of conformational changes of 
core spliceosomal components, stabilizing RNA-RNA- as well as 
RNA-protein interactions [reviewed in 51]. The only known 



enzymatic function of the NTC is a conserved U-box domain with 
E3-ubiquitin ligase activity of the Prp 1 9-protein [52,53]. It was 
conclusively shown that individual core components of the NTC, 
as Prp 19, Cefl/Cdc5 in S. cerevisiae and S. pombe, or the S. pombe 
SPF27 homologue Cwf7 independently affect splicing [42,50,54— 
59]. However, it is mostly unresolved to which extend single 
components of the NTC influence splicing of individual nitrons. 
We now show that one of the core components of the NTC, the 
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Figure 7. Deletion of the gene for a Numl -homologue in the 
ascomycete Aspergillus nidulans leads to a phenotype reminis- 
cent to the Jnuml phenotype in U. maydis. (A) The morphology of 
vacuoles is affected by the /4W4244-deletion. Spore suspensions of a 
wild-type (TN02A3) and the AN4244 deletion mutant were incubated 
for 18 hours at 37°C and the germinated filaments were CellTracker 
Blue (CMAC)-stained to visualize the vacuoles. Similar to the numl- 
deletion in U. maydis, the /4N4244-deletion in A. nidulans leads to an 
increased number of vacuoles. (B) Deletion of AN4244 affects nuclear 
positioning. Spore suspensions of A. nidulans wild-type (TN02A3) and 
the AN4244 deletion mutant were incubated for 12 and 18 hours, 
respectively, at 37"C, and the germinated filaments were DAPI-stained 
to visualize nuclei. In wild-type filaments, nuclei are evenly spaced, 
whereas in the mutant nuclei are often found in clusters. (C) The 
localization of Dynein (NudA) is dependent on AN4244. Spore 
suspensions of strains LZ12, expressing a GFP:NudA fusion protein, 
and LZ12zf/4N4244 were incubated for 18 hours at 37 C C. In LZ12 the 
GFP:NudA signal strongly accumulated in tips of germinating hyphae, 
whereas in mutant filaments, this apical accumulation was reduced; 
instead a stronger signal was observed in the cytoplasm. Scale bars: 
10 (im. 

doi:10.1371/journal.pgen.1004046.g007 
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SPF27 homologue Numl, has a global effect on splicing in U. 
maydis, as retention rates of nearly 80% of all introns were 
significantly increased in ram/-deletion mutants. We observed 
variations in intron retention between different transcripts in 
response to the Anuml mutation. Similar observations have been 
made in S. cerevisiae on a genome-wide scale for depletion of 
individual components of the spliceosome, including the NTC- 
components Prpl9 and Prpl7. The resulting splicing defects 
displayed considerable differences in different mRNAs, indicating 
that individual splicing factors affect a discrete set of genes [60] . It 
is well possible that structural alterations of the spliceosome 
complex after deletion oinuml result in altered recognition/affinity 
of specific introns, probably caused by their different GC content 
or secondary structures. Following the conclusions made in S. 
cerevisiae, we favor the possibility that deletion of numl leads to an 
altered integrity of the NTC, resulting in a reduced activity and/ or 
specificity and hence a diminished splicing efficiency. However, we 
cannot rule out that Numl may recruit splicing enhancer- or 
inhibitor-proteins that might affect splicing of only a subset of 
introns. 

Splicing defects can account for pleiotropic phenotypes 
in num 7-deletion mutants 

In addition to disturbed polarity we frequently observed 
delocalized septa in nam/-deletion strains; in addition, chitin- 
staining indicate structural differences to wild-type septa. Forma- 
tion of septa requires a highly coordinated and dynamic 
positioning of septins into ring-like structures [61], a process 
tightly regulated by the small GTPase Cdc42, its guanine 
exchange factor Donl and the protein kinase Don3 [62-64] as 
well as by formins or Cdc4, an essential light chain of the type II 
myosin motor protein [65]. The genes encoding the septins Cdc3 
(uml0503) and Sep3 (um03449), but also Cdc4 {umll848) display 
significandy increased intron retention rates (Table SI) and the 
expression of don3 (um05543) is reduced more than 2-fold in the 
nami-deletion strain (Table S4), which altogether might rationalize 
the formation of delocalized septa. 

We found that Numl affects splicing of several genes of the 
functional class "cell cycle", including genes for the catalytic and 
regulatory subunits of cyclin-dependent kinases, for mitotic cyclins 
as well as for y-Tubulin (Table S2). Similarly, several genes 
grouped in the functional class "DNA-repair" (e.g. for the DNA 
repair recombinase Rad51, the nucleotide excision machinery 
component Radl and the repair exonuclease Reel) show 
increased intron retention rates (Table S2). The synergistic effect 
of several genes with reduced splicing efficiency (and the resulting 
reduced activity) within the same functional category may be 
sufficient to explain the effects on the corresponding cellular 
processes. However, various components of the NTC have been 
found to be directiy involved in DNA damage repair. Initially, the 
resistance to the interstrand cross-linking reagent psoralen in a S. 
cerevisiae Prpl9 mutant strain [32] has been attributed to the 
splicing defect of the intron containing RAD 14 gene [66]. 
Subsequently, it has been conclusively shown that the DNA 
damage phenotype is independent from the splicing defect [42] . In 
addition, for Prpl9 several interaction partners have been 
identified which belong into the functional classes "cell cycle", 
"chromatin structure" and "DNA repair" rather than pre-mRNA 
processing. For example, the human Prpl9/Pso4 directly interacts 
with terminal deoxynucleotidyl-transferase (TdT), which is 
involved in DSB repair [67]. The human Cefl -homologue 
CDC5L interacts with the Werner syndrome protein WRN, a 
DNA helicase that functions during homologous recombination in 



response to DNA-damage [35] and also with the non-homologous 
end-joining (NHEJ) factor DNA PKcs [68]. 

Thus, it is well possible that also in U. maydis the increased UV- 
sensitivity of Anuml mutants results from a direct function of 
Numl in a DNA damage repair pathway. Similarly, disturbed cell 
cycle and cell division can be explained as indirect effect of a 
disturbed DNA-damage repair. Interaction of the human CDC5L 
with the cell cycle checkpoint kinase ATR ("ataxia-telangienctasia 
and Rad3-related") was shown to be required for activation of the 
ATR kinase as well as for further downstream effectors and 
mediators of ATR function as the checkpoint kinase CHK1 [69]. 
In U. maydis, the homologues of ATR 1 and CHK1 are known to 
be critical cell cycle regulators in addition to their decisive roles in 
the DNA-damage response [70,71]. 

We have shown previously that cell cycle and polarized growth 
in U. maydis is regulated by a cascade of transcription factors 
controlled by the bE/bW-heterodimer [2]. Several of the 
transcription factors of this b-cascade were found to have 
increased intron retention rates or a differential expression pattern 
in AB31 Anuml and might collectively contribute to the effects on 
cell cycle control and polar hyphal growth. For example, splicing 
efficiency of the hdpl gene, which encodes a homeodomain 
transcription factor required for filamentous growth and cell cycle 
control (Pothiratana and Kamper, unpublished data), was reduced 
from 83% in wild-type to 14% in the numi-deletion strain. More 
importandy, all four introns of rbfl are spliced less efficiently in the 
nami-deletion strain. As a consequence, protein levels of Rbfl, the 
master regulator for hyphal development, are significandy 
decreased in Anuml strains. In accordance with the reduced 
Rbfl protein levels in Anuml strains, we observed the down- 
regulation of rA/7-dependently expressed genes. Among these are 
the ri/7-dependendy expressed genes encoding the transcription 
factors Bizl and Hdp2, which were previously described as cell 
cycle regulators and factors for filamentous growth and pathogenic 
development [43; Scherer and Kamper, unpublished data]. 
Conclusively, we can assume that the reduced splicing efficiency 
of the rbfl gene and in turn the deregulation of rA/7-dependendy 
expressed genes contributes to hyphal morphology and the 
deregulated cell cycle in nami-deletion mutants. A direct proof 
of this hypothesis by replacing the rbfl open reading frame with a 
cDNA copy failed. The intron-free copy, although expressed to the 
same level as the intron-containing copy, was not able to provide 
wild-type Rbfl function, even when introduced into wild-type 
cells. Such a reduced expression of cDNA copies in relation to 
intron-containing genes has been described frequendy before, and 
in various cases this has been related to direct effects of introns to 
translational efficiency [reviewed in 72,73]. 

Numl has a specific role in cytoplasmic trafficking 
processes 

Aberrant or incomplete splicing events of the genes encoding 
motor proteins by themselves would be an explanation for the 
impaired endosomal trafficking processes observed in numl- 
deletion strains. However, none of the genes for the major 
transport-mediating motor proteins (kinl, kin3, dynl/2) contains 
introns or show altered expression levels in the A numl -mutant. 
Although we cannot rule out an effect caused by aberrant splicing 
of genes for unknown components of the transport machinery, we 
favor a specific role of the Numl protein in intracellular trafficking 
processes, based on the observation that (a) Numl is localized in 
distinct foci in the cytoplasm, implicating functions beyond its role 
as a splicing factor, and (b) Numl interacts with the Kinl motor 
protein, which is involved in long-distant transport processes 
during filamentous growth [9,12,14]. The physical interaction of 
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Kinl and Numl is in line with the similar phenotypes of the two 
mutant strains with respect to cell morphology, determination of 
polarity and septation, and morphology of vacuoles. Amml and 
likewise Akinl hyphae display impaired trafficking of early 
endosomes and reduced apical localization of Dynein, two 
interconnected phenotypes, as the transport of early endosomes 
relies on Dynein and the anterior transport of Dynein depends on 
Kinesin 1. In addition, the distribution of the direct Kinl cargo 
Mcsl is altered in both Akinl and Anuml mutants, which argues 
for a direct connection of Kin 1 and Num 1 . However, we can rule 
out that Numl is required for all Kinl -related processes, as (a) 
both mutations appear not to be epistatic, and (b) the numl and 
kinl mutant strains have overlapping, but not identical pheno- 
types; for example, haploid Akinl strains are defective in mating 
[20], while Anuml strains are not. This is anticipated, as each of 
the proteins participates in different cellular processes. The 
conventional kinesin is best known for its function in cellular 
transport processes. However, conventional kinesin has also been 
implicated in cell signaling, or to serve as bona fide, static anchor for 
mRNA-protein complexes or vesicular compartments, or to 
participate in the phosphorylation of cargo proteins to which it 
binds [reviewed in 74]. As outlined above, also SPF27 (or the 
NTC) participates in a wide spectrum of processes that, however, 
show no overlap with Kin 1 -dependent processes. 

Interaction of SPF27 homologues with conventional kinesin 
motor proteins is not specific to U. maydis, but appears to be 
conserved between basidio- and ascomycetes, as in Aspergillus 
nidulans the phenotype of deletion mutants of the SPF27 
homologue AN4244 or the Kinl homologue Kin A [48] are 
comparable to those observed for U. maydis. 

The association of a conventional kinesin motor protein with a 
component of the spliceosome-associated NTC has not been 
described so far and implies a novel function for SPF27- 
homologous proteins. The question is whether the Numl protein 
fulfills dual functions, (a) during splicing as an NTC-component in 
the nucleus and (b) during Kin 1 -dependent intracellular trafficking 
processes in the cytoplasm, or whether these two disparate 
mechanisms are connected by one single function. Discrete 
cytoplasmic functions may serve to regulate the activity of Kin 1 , 
as described in neuronal development for the interaction of 
Kinesin- 1 andJIP3, an interactor of the c-Jun N-terminal kinase 
[75]. Alternatively, Numl could function as a scaffold or adapter 
molecule, providing stability to the motor protein/cargo-complex 
and/or to the NTC. In both scenarios, deletion oinuml could alter 
the specificity /activity of Kin 1 -mediated transport processes, 
resulting in reduced motility of early endosomes and by that 
influence polar growth. We have not been able to visualize Numl 
proteins moving or co-migrating with Kinl along microtubule 
tracks. As GFP-Kinl shows a strong cytoplasmic background, 
visualization of Kin 1 molecules moving on microtubules has not 
been achieved yet [14,76]. Taken into account that the interaction 
of Kin 1 and Num 1 appears to be transient, we cannot rule out that 
Kinl and Numl may interact on microtubule tracks. 

A possible process that could couple both nuclear and 
cytoplasmic functions of Numl is the cytoplasmic transport of 
mRNAs. In U. maydis, the RNA binding protein Rrm4 was found 
to be instrumental for bidirectional transport of mRNAs in 
filaments [47]. Intriguingly, cytoplasmic transport of Rrm4 is 
dependent on Kin 1 , Kin3 and Dynein, and recently it was shown 
that Rrm4 co-localizes with early endosomes [7,77]. There is a 
remarkable overlap of the phenotypes of rrm4-, numl- as well as 
fani-deletion mutants with respect to polarized growth [14,78], 
which raises the possibility of a concerted function of all three 
proteins within the same process. Indeed, the motility of early 



endosomes associated with Rrm4:eGFP fusion proteins is drasti- 
cally reduced in Afe'ni -mutants [7,14] and, as we show here, also 
in rami-deletion strains. 

In Drosophila oocytes, it has been shown that components of the 
exon junction complex (EJC) are essential for transport of oskar 
transcripts [79], The EJC marks exon/exon-transitions during 
splicing, remains bound to the processed mRNA and serves as a 
quality control mechanism for the splicing reaction during 
translation [reviewed in 80]. Reminiscent to mRNA transport in 
U. maydis, transport of oskar-mRNA is also dependent on the 
conventional Kinesin 1 [81]. 

Unlike the components of the EJC that remain bound to mRNA 
during nuclear export, the NTC is thought to disassemble after the 
splicing reaction. However, recendy it has been shown in S. 
cerevisiae that the NTC is required for efficient recruitment of 
TREX to the transcriptional machinery. This is mediated by the 
C-terminus of the NTC-component Syflp, which connects TREX 
to Rpbl, the largest subunit of RNA-polymerase II [82]. TREX is 
a conserved protein complex ("transcription export"), which 
facilitates both transcription elongation as well as mRNA-export 
processes [83]. Yralp, a component of the TREX complex, is an 
RNA-binding protein and interacts with the Mex67/Mtr2- 
receptor complex, which in turn act together with components 
of the nuclear pore complex (NPC) to guide the transcripts into the 
cytoplasm [84-86]. Thus, it appears well possible that NTC 
components remain tethered to the spliced mRNA. In this context 
it is interesting to note the recently discovered association of 
Kinesin 1- and Dynein to fungal nuclear pores, which mediate 
NPC-motility and are thus required for chromosome organization 
and nucleo-cytoplasmic transport [76]. This result corroborates a 
functional coupling between the nuclear transcription machinery 
and cytoplasmic transport processes in U. maydis. 

The major finding of our study is that a component of the 
conserved spliceosome-associated NTC is involved in cytoplasmic 
trafficking processes. Due to the interaction of Numl with the 
Kinl and the functional connection to the RNA-binding protein 
Rrm4, it is tempting to speculate that Numl could have a function 
during intracellular mRNA-transport processes. We hypothesize 
that one of the functions of the Num 1 -protein is to participate in 
the coordination of pre-mRNA splicing and the formation and 
export of mRNP-particles out of the nucleus. The transcripts 
would be transferred into the cytoplasm where they could be 
passed on to the microtubule-based transport machinery via the 
interaction of Numl with Kinl motor proteins that were found to 
be associated with the nuclear pore complex (Figure 8). In this 
model, additional adapter-proteins that act as connecting proteins 
between the nuclear splicing machinery and the cytoplasmic NPC- 
bound motor proteins cannot be excluded and remain to be 
identified. 

Materials and Methods 

Strains and growth conditions 

Escherichia coli strain TOP 10 (Invitrogen) was used for cloning 
purposes. Growth conditions and media for the cultivation of E. 
coli followed the protocols described previously [87]. The 
Saccharomyces cerevisiae strain AH 109 (Clontech) was used for yeast 
two-hybrid interaction studies. S. cerevisiae cells were grown in 
YPDA complete medium, or on minimal medium (SD) supple- 
mented with the dropout-mix needed for selection, as described in 
the Clontech Matchmaker™ GAL4 Two-Hybrid System 3 
Manual (http://www.clontech.com). Ustilago maydis cells were 
grown in YEPSL [88], CM (complete medium) supplemented 
with 1% glucose (CM-G) and 1% arabinose (CM-A), respectively, 
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Figure 8. Hypothetical model for function of Numl coordinat- 
ing pre-mRNA splicing with nuclear pore complex-dependent 
export of mRNP-particles and microtubule-based mRNA-trans- 
port. Num1 is a component of the spliceosome-associated Prp19/ 
CDC5L-complex (NTC) and has a global effect on splicing efficiency. 
NTC-components interact with the TREX ("transcription-coupled 
export")-complex [82], which mediates the export of mature transcripts 
into the cytoplasm via the nuclear pore complex (NPC) [83-86], Adapter 
proteins connect the mRNA to NPC-bound export receptors. The 
Kinesin-1 motor protein has been shown to be associated to the NPC 
and promotes nucleo-cytoplasmic transport processes [76]. mRNAs 
within the hyphae are distributed via the RNA-binding protein Rrm4, 
which shuttles associated with early endosomes along the polar 
microtubule cytoskeleton [7,77]. Motility of the Rrm4 associated early 
endosomes depends both on Kinl [7,14] and Num1. Proteins labeled 
with (?) depict yet unidentified adapter-proteins that might function as 
connectors between the nuclear splicing machinery and the cytoplas- 
mic NPC-bound Kinl. As Num1 is a component of the NTC complex and 
interacts with Kinl, it might fulfill a function during mRNA export 
processes, possibly by connecting exported mRNP to the cytoplasmic 
transport machinery via its interaction with Kin 1 . 
doi:1 0.1 371 /journal.pgen.1 004046.g008 

or NM (nitrate minimal medium) [89] at 28°C. Solid media 
contained 2% agar. The induction of hyphal growth in AB31- and 
AB33 derivatives was done as previously described in [24]. Plate 
mating assays followed the protocol of [24]. Mating assays in liquid 
culture for the generation of dikaryotic hyphae were carried out as 
described in [90]. Pheromone stimulation of U. maydis cells was 
performed according to the protocol of [9 1] . Aspergillus nidulans was 
grown in CM (complex medium) [92] or MM (minimal medium) 
[93]. 

All U. maydis and A. nidulans strains used in this study are listed in 
Tables S7 and S8, respectively. 

Plant infections were performed as described in [94]. For 
infection studies with U. maydis, the maize cultivar Early Golden 
Bantam was used under controlled conditions in a CLF Plant 
Climatics GroBank with a 16 h (28°C)/8 h (22°C) day/night 
rhythm. Disease symptoms were evaluated according to the 
disease rating criteria published previously [23]. 



Yeast two-hybrid analyses 

Screening for Num 1 -interacting proteins was performed as 
previously described [95], using the Matchmaker III system 
(Clontech). Plasmid pGBKT7-Numl was generated by PCR 
amplification of the numl open reading frame (ORF) from 
genomic DNA (strain 521 [23]), introducing two incompatible 
SfH restriction sites at the 3' and 5' ends and subsequent ligation 
into the respective sites of pBGKT7 (Clontech). Plasmids pGAD- 
Prpl9 and pGAD-Cefl containing full-length as well as N- and C- 
terminal variants were generated by PCR amplification using 
either genomic DNA (for cefl) or a full-length cDNA library (for 
prpl9) and subsequent ligation into pGAD-DS (Dualsystems 
Biotech) via 5/3 sites. Oligonucleotide sequences are given in 
Table S9. 

DNA and RNA procedures 

Molecular techniques followed the protocols described in [87]. 
DNA from U. maydis was isolated according to the protocol given 
in [96]. Transformation procedures were performed as described 
in [97]. DNA from A. nidulans was isolated according to [98] and 
transformations were carried out as described in [99]. 

For gene deletions in U. maydis, a PCR-based approach [100] 
was used. Open reading frames were replaced by either 
hygromycin, carboxin or nourseothricin resistance cassettes. 
Similarly, for gene deletions in A. nidulans, a pyridoxine 
auxotrophic marker was used. C-terminal eGFP, RFP or 3xHA 
gene fusion constructs were generated in U. maydis via homologous 
recombination following the protocol given in [101], using SJil 
cassettes from plasmids pUMA647 (eGFP) [7], pUMA738 [47] 
(RFP) or pUMa791 or pUMA792 (Feldbrilgge, unpublished) for 
3 xHA fusions. Fusion constructs were subcloned in pCRII TOPO 
(Invitrogen), PCR-generated linear DNA was used for transfor- 
mation of U. maydis. Homologous integration of the constructs was 
verified by Southern analysis. All sequences of oligonucleotides 
used for PCR are listed in Table S9. 

Total RNA was extracted using Trizol reagent (Invitrogen) 
according to the manufacturer's instructions. RNA samples to be 
used for real-time qRT-PCR or mRNA-Seq analysis were column 
purified (RNeasy; Qiagen) and the quality was checked using a 
Bioanalyzer with an RNA 6000 Nano LabChip kit (Agilent). 
Quantitative RT-PCR analysis for pcna, rpb3, ubil, rho3 and rbfl 
was conducted with oligonucleotide pairs that specifically detected 
spliced and unspliced RNA species, respectively: forward primers 
used for quantification of spliced transcript were specific for exon/ 
intron borders, whereas forward primers for the unspliced 
transcripts were placed within the intron sequence. Ratios of 
spliced vs. unspliced RNA of AB31 (wild-type) samples were set to 
1.0 and compared to the ratios obtained with the AB31zl™ml. 
qRT-PCR experiments were done in three independent biological 
and two technical replicates each and followed the protocols given 
in [102] and [103]. 

DNA content was measured by flow cytometry. FACS analyses 
were performed according to [104] with the exception that 10 5 
cells were harvested for each time point and analyzed with a BD 
Biosciences FACSVerse flow cytometer using the FACSDiva 6.2 
software. For each acquisition, 2*10 events were measured at an 
average flow rate of 100-500 events per second. 

mRNA-sequencing 

mRNA was sequenced starting from total RNA isolated from 
strains AB31 and AB'HAnuml on an Illumina HiSeq 1000 using 
the TrueSeq mRNA sequencing kit according to the manufactur- 
er's instructions. Three independent libraries for each of the two 
strains were generated (Olivier Armant, Insitute for Toxology and 
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Genetics, KIT Campus North), using three independently grown 
cultures (biological replicates). Each library yielded an average of 
60 million paired-end reads of 2x56 bp, which equals about 
20xl0 9 bp sequences for each strain. Raw paired-end reads 
(2x56 bp) were aligned with TopHat 1.3.2 [105] to the Ustilago 
genome assembly provided by the MIPS Ustilago maydis Database 
(MUMDB, http:/ /mips. helmholtz-muenchen.de/genre/proj/ 
ustilago). TopHat was supplied with a gtf file derived from 
MUMDB's p3_t237631_Ust_maydi_201 10629.gif, but not re- 
stricted to known transcripts. Aligned reads were counted using 
custom Python scripts using the pysam library (http://code. 
google. com/p/pysam/). Gene expression values were derived by 
counting read pairs that were completely inside the exon regions of 
MUMDB-genes and differential expression was assessed with 
DESeq 1.6.1 [106] and BaySeq 1.8.1 [107] and filtered for an 
adjusted p-value^0.05 and a more than two fold change. 

For comparison and visualization of splicing efficiency, read 
pairs overlapping the respective intron/exon with their first 
segment by at least three base pairs were counted and normalized 
to fragments per kb per million (FPKM). Intron retention was 
determined by calculating the ratio of introns FPKM to the gene's 
overall exon FPKM. Exon FPKM included all exons with more 
than 10 bp. Only genes with an FPKM above 10 were considered. 
A small number of introns with a retention rate above 0.9 were 
excluded, based on the assumption that they base on false gene 
models. Sequencing data was deposited at the EBI ArrayExpress 
Database (E-MTAB-1300). 

Co-immunoprecipitation analysis 

For in vitro protein expression and co-immunoprecipitation of 
Numl and Kinl, the TNT quick coupled transcription/transla- 
tion system (Promega) was used according to the manufacturer's 
protocol. For expression of Myc-tagged Num 1 protein, plasmid 
pGBKT7-Numl was used and HA-tagged Kinl protein was 
generated using the pGAD-Kinl 65 o_9 68 derivative identified in the 
yeast two-hybrid screen. All subsequent steps were carried out as 
described previously [103]. 

For in vivo co-immunoprecipitation, U. maydis sporidia were 
grown in 100 ml GM-G to an OD 600 of 0.8. Cells were washed 
once with PBS buffer, resuspended in 1 ml PBS supplemented 
with Complete proteinase inhibitor cocktail (Roche), frozen in 
liquid nitrogen and homogenized in a Retsch mill for 10 minutes 
at 30 Hz. The cell lysate was incubated with 40 |J,g monoclonal 
anti-HA coupled agarose (Sigma-Aldrich) on a rotating wheel at 
4°C over night. Agarose beads were washed three times in PBS 
prior to resuspension in 20 u.1 Laemmli buffer. Samples were 
boiled for three minutes and separated by SDS-PAGE. For co- 
immunoprecipitation of Numl and Kinl, a crosslinking reaction 
with dithiobis[succinimidylpropionate] (Thermo Fisher) was con- 
ducted according to [108] using 300 ml of cells grown in CM-G to 
an OD 600 of 0.8. 

Proteins were transferred to PVDF nitrocellulose membranes in 
a semidry blot chamber. Western blots were probed with 
monoclonal anti-HA, anti-c-Myc, and anti-GFP (Sigma-Aldrich) 
antibodies. Horseradish peroxidase-conjugated anti-mouse or 
anti-rabbit IgG (Promega) was used as secondary antibody, and 
an ECL system was used for protein detection. 

Microscopy, image processing and quantitative analysis 

For microscopic analyses, logarithmically growing U. maydis cells 
were taken from liquid cultures grown in CM-G medium. For the 
induction of hyphal growth, cells were shifted to CM-A or NM to 
induce the P c , g r or P narl - promoters, respectively, for 12-14 hours. 
For microscopy of A. nidulans germlings and young hyphae, MM 



on cover slips was inoculated with a small amount of spores and 
incubated for 12-18 hours at 37°C. Cells were then placed on top 
of a 2% agarose cushion placed on the microscope slide and 
immediately observed using an Axioimager Zl microscope 
equipped with an Axiocam MRm camera (Carl Zeiss). Standard 
filter sets for DAPI, GFP, CFP and Rhodamine were used for 
epifluorescence analysis. 

Nuclei were stained with DAPI Vectashield H-1200 (Vector 
Laboratories); chitin was stained with 2 Hg/ml Calcofluor White 
(Sigma-Aldrich) or with Congo Red solution (1 |jl/ml; stock: 
1 mg/ml in H 2 0). Congo Red was incubated at 25°C for 
10 minutes on a turning wheel and washed twice with CM prior to 
microscopic analysis. Fungal cell walls in general were stained with 
1 mg/ml WGA/Fluorescein (Invitrogen) in PBS. For the visual- 
ization of vacuoles, growing cells were incubated in CM 
supplemented with 10 u,g/ml CellTracker Blue (7-amino-4-chlor- 
omethyl-coumarin, CMAC) (Invitrogen) for 30 minutes at 28°C, 
washed twice with CM and then subjected to microscopy. 
Chlorazole Black E staining was performed according to [24] . 

Endosome motility was measured in image sequences of 50 
frames, taken with an exposure time of 500 msec. Only organelles 
that were moving for a distance of at least 5 |0.m were considered. 
The resulting movies were converted into kymographs using 
ImageJ software [109]. For quantification of Mcsl signal intensity, 
the corrected total cell fluorescence (CTCF) of mid-size growing 
buds was calculated with ImageJ software (CTCF = integrated 
density - (area of selected cellxfluorescence of background 
reading)). 

For quantification of the microtubule orientation, cells were 
grown to mid-log phase in YEPSL medium. Cells were analyzed 
via 45 sec time-lapse recordings (four frames each 15 sec). Pictures 
were taken with an exposure time of 90 msec for Tub 1 -GFP and 
300 msec for PebhRFP signals. Movement of Pebl:RFP signals 
was scored at a screen as described in [46] . All image processing, 
including adjustment of brightness, contrast and gamma-values 
was performed with the AxioVision and ZEN software (Carl 
Zeiss), respectively. 

Phylogenetic analysis 

For comparative phylogenetic analysis of Num 1 , 65 sequences 
with the highest similarity to Numl were obtained by BLASTP 
analysis. The sequence of Saccharomyces cerevisiae Snt309p was 
included as outgroup. Sequences were aligned with MAFFT 
version 6 [110] using the global alignment G-INS-i. A phyloge- 
netic tree was calculated using the minimum linkage clustering 
method (http:/ / align.bmr.kyushu-u.ac.jp/ mafft/ online/ server/). 

Fig Tree 1.3 (http://tree.bio.ed.ac.uk/software/figtree/) was 
used to visualize the Nexus formats of the MAFFT results. 

Accession numbers 

Sequence data from this article can be found at the Munich 
Information Center for Protein Sequences (MIPS) Ustilago maydis 
database (http://mips. helmholtz-muenchen.de/genre/proj/ 
ustilago/), the Aspergillus genome database (http://www.aspgd. 
org/) and the National Center for Biotechnology Information 
(NCBI) database under the following accession numbers: numl 
\um()1682), XP_757829; pip 19 (uml0027), XP_756259; cefl 
\um04411), XP_760558; rbfl {um()3172), XP_759319; pcna 
\um05403), XP_761550; rpb.3 (um03550), XP_759697; ubil 
{um02440), XP_758587; rho3 (um04070), XP_760217; rad51 
\um03290), XP_759437; reel (umll016), XP_759527.1; mms2 
\uml0097), XP_756467.1; cksl (um()3210), XP_759357; cdkl 
[uml0705), AAP94021.1; hdpl (uml2024), XP_761909.1; hdp2 
{um04928), XP_761075; bid (um02549), XP_758696; um02704 
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(related to allantoate permease), XP_758851; uml2105 (probable 
PUP3 - 20S proteasome subunit), XP_756658.1; actin (umll232), 
XP_762364; eIF2b (um04869), XP_761016; kinl (um04218), 
XP_760356; dyn2 (um04372), XP_760519; yupl (um05406), 
XP_761553; tubl {um01221), XP_757368; AM4244, 
XP_66 1848.1; nudA [AN0118], XP_657722.1; kinA (AN5343), 
XP_662947.1.; the gene model for uml5049 (related to PUF3 
transcript specific regulator of mRNA degradation) is implement- 
ed only at MUMDB (http://mips.helmholtz-muenchen.de/ 
genre/proj/ustilago), as the predicted protein is composed of 
two fragments located on adjacent contig ends and, according to 
the sequence gap, part of the protein may be absent. 

Supporting Information 

Figure SI Alignment of Numl to SPF27 homologues. ClustalW 
was used to align proteins related to Numl from S. reilianum 
(Srl2752, NCBI accession number CBQ71896), U. hordei 
(UHOR_02497, CCF53376) C. cinerea (Numl, BAC78620), A. 
nidulans (AN4244, XP_66 1848.1), A. thaliana (Mos4, 
AT3G18165.1), mouse (SPF27, BAB31409) and human 
(BCAS2/SPF27, CAG46834). Identical or similar residues are 
highlighted in black and grey, respectively. The nuclear localiza- 
tion signal (NLS) in the U. maydis sequence is highlighted in green. 
The red box indicates the conserved BCAS2-domain. 
(PDF) 

Figure S2 Phylogeny of Numl homologues. For comparative 
phylogenetic analysis of Numl, 65 sequences with the highest 
similarity to Numl were obtained by BLASTP analysis. The 
sequence of Saccharomyces cerevisiae Snt309p was included as out- 
group. Sequences were aligned with MAFFT version 6 [110] using 
the global alignment G-INS-i. A phylogenetic tree was calculated 
using the minimum linkage clustering method. 
(PDF) 

Figure S3 Deletion of numl does not affect growth in axenic 
culture, mating and dikaryon formation or proliferation in planta. 
(A) Growth curves of AB31 and AB'ilAnuml deletion strain in 
YEPSL complex- (black) and glutamin/glucose-containing mini- 
mal medium (grey). The numl -deletion has no influence on growth 
in axenic culture. (B) SG200 (control) and SG200Anuml sporidia 
were grown in glucose-containing CM-medium. No obvious 
phenotype could be observed, indicating that the A numl '-mutant 
is not impaired during growth in axenic culture. (C) Mating assays 
on charcoal-containing CM-glucose medium by co-spotting of 
strains indicated. FBI (albl) and FB2 (a2b2) strains are included as 
wild-type controls. Formation of dikaryotic aerial hyphae is visible 
as white mycelium. The ram/-deletion has no influence on the 
formation of the dikaryon. (D) Pathogenicity of individual numl- 
deletion strains. Disease rating of maize seedlings seven days post 
inoculation with U. maydis strains SG200, SG200/f««mi, FBI xFB2 
(wild-type crosses) and FBI Anuml xFB2zfnam7-derivatives. #1, 
#2 indicate independently obtained deletion mutants. Bars 
represent the percentage of infected plants with the symptoms 
indicated in the legend. N corresponds to the total number of 
plants infected. (E) Maize plants were inoculated with SG200 and 
SG200Anuml . Infected leaves were stained with Chlorazole Black 
E. Samples were taken three and five days post inoculation, 
respectively. Clamp-like structures (arrowheads) could be observed 
in both strains, indicating fungal proliferation within the plant 
cells. Scale bars: 10 nm. 
(PDF) 

Figure S4 Phenotype of conjugation hyphae of nami-deletion 
strains.. Strains FBI and FBlAnuml were treated with synthetic 



a2-pheromone (2,5 u,g/ml) for six hours in glucose-containing 
CM-medium. DMSO was used as a solvent for the a2-pheromone 
and served as negative control. Both FBI and FBlAnuml respond 
to the pheromone with the formation of conjugation tubes. Note 
the characteristic A numl -mutant phenotype, such as branches and 
irregular hyphal growth. Scale bars: 10 |0.m. 
(PDF) 

Figure S5 Phenotype of dikaryotic hyphae of numl -deletion 
strains. Compatible FBI- and FB2-strains were used to monitor 
the zlnami-phenotype in dikaryotic hyphae. To demonstrate 
successful fusion events, FBI -strains express mitochondrial ma- 
trix-targeted RFP (mtRFP), and FB2-strains express mitochondrial 
matrix-targeted GFP (mtGFP), respectively, under control of the 
inducible P c „,;-promoter. Fusion of the strains was induced in 
arabinose-containing liquid charcoal medium overnight. Depicted 
are DIC-, GFP- and RFP-signals. Merged images show both 
fluorescent signals within the same hypha, indicating fusion of two 
sporidia in compatible FBI xFB2 wild-type crossings (upper panel) 
as well as in compatible FBlAnuml xFB2Anuml crossings (lower 
panel). The two initial cells are marked with an asterisk. Note the 
characteristic Anuml -rmiVAXiX phenotype (thicker, irregular hyphal 
growth and bipolar growth of the sporidia). Scale bars: 10 urn. 
(PDF) 

Figure S6 The zframi-deletion causes chitin accumulations and 
the insertion of unusual septa. Induction of the bEl/bW2 
heterodimer in wild-type cells (AB31) induces a G2 cell cycle 
arrest, resulting in hyphae with a single nucleus positioned in the 
tip compartment. 18% of AB3lAnuml mutant hyphae contain 
more than one nucleus within the cell. Filaments were stained with 
DAPI to visualize nuclei and Fluorescein-conjugated wheat germ 
agglutinin (WGA), which specifically binds to jV-acetylglucosamine 
[27]. In AB31 wild-type cells, chitin was intensively stained in the 
growth region within the hyphal apex as well as in the basal 
retraction septum (marked with arrowhead). In ABSlAnuml 
mutants, however, chitin accumulations along the hyphae (marked 
by asterisks) as well as septa-like chitin rings were observed 
frequently (middle and lower panel). In many cases, these septa were 
not visible in the DIC channel, indicating that these structures may 
not represent true septal cell walls. Scale bars: 10 [Lm. Septa in 
AB31 Anuml were also visualized by Calcofluor (Figures IE and 6B) 
and Congo Red (Figure 3A). In contrast to septa in the dikaryotic 
hyphae, septa in AB3lA?iuml were not always visible in the DIC 
channel, indicating that these septa or septa-like structures may be 
structurally different from wild-type septa. 
(PDF) 

Figure S7 Complementation analyses of nami-deletion strains. 
Filaments of AB31 and AB3lAnuml, expressing NumlxGFP 
fusion proteins under control of the native P num r or the inducible 
P „ ; -promoter, were analyzed 12-14 hours after induction of the 
bEl/bW2-heterodimer and, when present, the NumlxGFP 
fusion proteins. (A) In AB31, NumlxGFP fusion proteins localize 
predominantly to the nucleus. (B) Complementation of the numl- 
deletion with a NumlxGFP fusion protein restores the wild-type 
phenotype and the protein localizes predominantly to the nucleus. 
(C) A Numl:3eGFP fusion protein with mutated nuclear 
localization signal (NLS, shown in D) does not complement the 
Anuml phenotype; filaments show the typical curved and bipolar 
morphology of the nami-deletion. In addition to the GFP-signal in 
the nucleus, a strong signal is observed in the cytoplasm. Scale 
bars: 10 urn. (D) Schematic representation of the Numl protein 
depicting the basic amino acids within the NLS that were replaced 
by alanine-residues. 
(PDF) 
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Figure S8 Prpl9- and Cefl -fusion proteins are functional and 
are localized in the nucleus. prpl9 and cefl were expressed as C- 
terminal RFP- or HA-tagged versions under control of their 
endogenous promoters in strain AB31 to ensure native expression 
levels. Hyphal growth was analyzed 1 4 hours after induction of the 
bEl/bW2-heterodimer. (A) Shown are DIG- and RFP-fluores- 
cence signals of strains UNK208 (AB31 numl:3egfp:hyg R , 
prpl9:rfp:nai i ) and (B) UMO10 (AB31 numl:3egfp:hyg R , cefl:rfp:naf). 
As prpl9 and cefl are both essential genes (see Supporting Text SI), 
the viability of the strains proofs the functionality of the fusion 
proteins. In addition, no altered phenotypes with respect to 
filamentous growth were observed, corroborating that both 
Prpl9:RFP as well as Cefl:RFP fusion proteins are functional. 
Scale bars: 10 u,m. 
(PDF) 

Figure S9 Cell cycle analysis of AB31 vs. AB3\Anuml. (A) FACS 
analysis of AB31 and AB'ilAnuml after induction of the bEl/ 
bW2-heterodimer in arabinose-containing CM-medium. Samples 
were taken at time points indicated. The histograms show the 
DNA-content measured by FACS analysis. Relative fluorescence 
intensities are given on horizontal axes, vertical axes reflect cell 
numbers. In AB31, induction of the active bEl/bW2-combination 
resulted in an enrichment of cells with a 2C DNA-content, 
indicative of a G2 cell cycle arrest. In AB'il Anuml, a higher 
amount of cells in G 1 was observed, which implicates defects in 
cell cycle regulation. (B) Microscopic analysis of the filamentous 
growth of the strains upon b- induction. Scale bars: 10 |im. 
(PDF) 

Figure S10 Control of bW-, rbfl- and numl-gene expression in 
samples used for RNA-Seq-analysis. RNA was isolated from 
strains AB31 and AB'ilAnuml in three biological replicates each, 
eight hours after induction of the bEl/bW2-heterodimer in 
arabinose-containing minimal medium. Expression of bW-, rbfl- 
and numl was monitored by qRT-PCR analysis. Gene expression 
is shown relative to the highest expression value. Actin and eIF2b 
were used for normalization. Shown are the mean values two 
technical replicates. Error bars represent the SD. 
(PDF) 

Figure SI 1 Verification of mRNA-Seq using qRT-PCR. (A) Plot 
depicting splicing efficiency of genes based on the RNA-Seq 
analysis. Plotted are the FPKM values (fragments per kilobase of 
sequence per million fragments mapped) across the genomic 
region indicated (coordinates in nucleotides) of three independent 
RNA-Seq experiments for AB31 wild-type (blue lines) and 
AB'ilAnuml (red lines), respectively. Exons (E) and introns (I) 
are indicated. (B) Schematic representation of exon-intron 
structures. Oligonucleotides used for qRT-PCR are depicted as 
arrows and allow the discrimination against spliced (primer 
depicted in red) and unspliced transcripts (primer depicted in 
green). (C) Gene expression analyses to verify expression of spliced 
vs. unspliced transcripts of the indicated genes using qRT-PCR. 
RNA was isolated from AB31 and AB'ilAnuml eight hours after 
induction of the bE/bW-heterodimer. Actin and eIF2b were used 
for normalization. Depicted are the ratios of spliced vs. unspliced 
transcripts of three independent biological replicates. Mean values 
of two technical replicates each are shown. 
(PDF) 

Figure SI 2 Individual genes show different alterations in intron 
retention in AB'ilAnuml . Splicing efficiency of the indicated genes 
based on the RNA-Seq analysis. Depicted are examples for genes 
(A) where splicing is not affected, (B) the intron is not spliced at all, 
or (C) different introns show dissimilar retention rates. Plotted are 



the FPKM values (fragments per kilobase of sequence per million 
fragments mapped) across the genomic region indicated (coordi- 
nates in nucleotides) of three independent RNA-Seq experiments 
for AB31 wild-type (blue lines) and AB3\Anuml (red lines), 
respectively. Exons (E) and introns (I) are indicated. 
(PDF) 

Figure S13 Reduced splicing efficiency in AB31 Anuml in genes 
with function in DNA repair and cell cycle. Splicing efficiency 
based on RNA-Seq analysis. Depicted are examples for genes 
grouping in the functional category (A) "DNA repair" and (B) 
"cell cycle". Plotted are the FPKM values (fragments per kilobase 
of sequence per million fragments mapped) across the genomic 
region indicated (coordinates in nucleotides) of three independent 
RNA-Seq experiments for AB31 wild- type (blue lines) and 
AB'ilAnuml (red lines), respectively. Exons (E) and introns (I) 
are indicated. 
(PDF) 

Figure S14 Splicing defect of the rbfl-gene in ABilAnuml. (A) 
Schematic view of rbfl intron/ exon structure and primers used for 
RT-PCR (blue) and qRT-PCR (red, green and black). (B) RNA 
samples of AB31 and AB31 Anuml were isolated eight hours after 
induction of the bEl/bW2-heterodimer. RT-PCR analysis was 
performed on the rbfl open reading frame and PCR products were 
separated on a 2% TAE-agarose gel. In AB31 cDNA, the spliced 
version of rbfl (430 bp) is predominantly detected. In AB31 Anuml 
more bands spanning different sizes from 430 bp (all four introns 
spliced) to 1192 bp (all introns retained) are observed, validating 
the splicing defect of the rbfl-gene observed in the RNA-Seq 
analysis. (C) qRT-PCR analysis to investigate the expression of 
spliced vs. unspliced ri/7-transcripts. Primers depicted in (A) were 
used to distinguish spliced (primer depicted in red) from non- 
spliced (primer depicted in green) introns. Actin and eIF2b were 
used for normalization. Depicted are the ratios of spliced vs. 
unspliced transcripts of three independent biological replicates. 
Mean values of two technical replicates each are shown. 
(PDF) 

Figure S15 Splicing efficiency and expression profile of the 
Rbfl -target genes hdpl, hdp2, bizl is altered in AB31 Anuml. Shown 
are the expression profiles based on the RNA-Seq analysis. Plotted 
are the FPKM values (fragments per kilobase of sequence per 
million fragments mapped) across the genomic region indicated 
(coordinates in nucleotides) of three independent RNA-Seq 
experiments for AB31 wild-type (blue lines) and AB31 Anuml 
(red lines), respectively. Exons (E) and introns (I) are indicated. (A) 
Expression of the hdpl-gene is increased about two-fold, and 
splicing efficiency of both introns is reduced in ABilAnuml. (B) 
and (C) Expression levels of hdp2 as well as bizl are significantly 
reduced in ABilAnuml (hdp2: -3,1-fold, p = 5.28 * 10" 36 bizl: 
— 3,4-fold, p = 5.34 * 10 62 ; significance according to DEseq 
[106]). 
(PDF) 

Figure S16 Introns are essential for full function of rbfl. (A) and 
(B) An intron-free derivative of rbfl was expressed in AB31, 
replacing the endogenous rbfl-gene. cDNA was isolated from 
strains AB31 and UT062 (AB31 rbfl ::rbfl cDNA) 10 hours after 
induction of filamentous growth in arabinose-containing CM- 
medium. qRT-PCR analysis was performed to verify unaltered 
expression of bE and rbfl (A). The intron-free rA/7-derivative is not 
able to induce rA/7-target genes (B). Gene expression is shown 
relative to the highest expression value. Actin and eIF2b were used 
for normalization. Mean values of two technical replicates each 
are shown. (C) Schematic representation of the construct used in 
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this study. rbfl-cDNA was expressed in the genomic rbfl-locus 
under control of the endogenous P r 4yj-promoter. 500 bp of the 3'- 
UTR were included as terminator. A carboxin resistance-cassette 
was used for selection. (D) SG200-derivatives were spotted on 
charcoal-containing CM-medium and incubated at 22°C for 3 
days. The formation of filamentous hyphae is visible as white 
mycelium. In contrast to the complementation with a wild-type 
rbfl -copy [2], the respective intron-free rA/7-derivative, which 
replaced the endogenous rbfl, was not able to complement the 
/(r^/7-phenotype. Magnified images of the colony margins were 
taken with a binocular assisted camera. (E) Maize seedlings were 
infected with the indicated SG200 derivatives. Disease rating and 
tumor formation was monitored seven days after infection. Bars 
represent the percentage of infected plants with symptoms 
indicated in the legend. N corresponds to the number of plants 
infected. The intron-free rA/7-derivative leads to a reduction in 
tumor formation. (F) and (G) Microscopic analysis of the indicated 
strains 10 hours after induction of filamentous growth in 
arabinose-containing CM-medium. Phenotypes were grouped 
into four different categories. Whereas 95% of AB31 wild-type 
cells grew as long, straight filaments, the intron-free rbfl -derivative 
was not able to fully complement the zlrA/7-phenotype. Integration 
of the ri/7-cDNA construct into the AramZ-mutant background 
hence did not reduce the observed phenotypes of the numl- 
deletion strain. (G) Representative microscopic images of the 
indicated strains. Shown are DIC-images. Scale bars: 10 urn. 
(PDF) 

Figure S17 Characterization of the interaction domain of Numl 
for the Kinl motor protein. N- and C-terminally truncated 
fragments of the Numl -protein were tested for interaction with 
Kinl in the yeast two-hybrid system. Left panel: Schematic 
depiction of the constructs used in this study. Numbers represent 
the regions of the Numl -constructs in amino acids. The conserved 
BCAS2-domain is depicted in dark grey. After cloning into 
pGBKT7 (Clontech), the constructs were co-transformed with the 
pGAD-Kinl g5o_968 fragment isolated in the yeast two-hybrid 
screen. Serial 10-fold dilutions of cell suspensions of the respective 
transformants were spotted on selective media (SD -LW) and 
stringent media (SD -LWHA), indicating an interaction. Plates 
were incubated for two days at 30°C. The Numl interaction 
domain resides between amino acids 67 and 147 (highlighted in 
red). pGBKT7-p53 with pGADT7-T served as positive control, 
pGBKT7-Lam with pGADT7-T (Clontech) served as negative 
control. 
(PDF) 

Figure S18 The double deletion of Numl and Kinl has additive 
effects on pathogenicity as well as hyphal morphology. (A) Effects 
of the double deletion on polar hyphal growth in strain 
AB31 Anuml Akinl . Filaments of were analyzed 12 hours after 
induction of the bEl/bW2-heterodimer. Compared to AB31 Akinl 
and AB31 Anuml hyphae, which predominandy grow as bipolar 
filaments, the double deletion has additive effects on filament 
morphology, polarity and branching. (B) To visualize nuclear 
positioning and septation, filaments of AB31 Anuml Akinl were 
treated with DAPI- and WGA/FITC. Scale bars: 10 um. (C) 
Pathogenicity of individual numl- and £m/-deletion strains. Maize 
seedlings were infected with the indicated SG200 derivatives. #1, 
#2, #3 indicate independendy obtained deletion strains. Disease 
rating and tumor formation was monitored seven days after 
inoculation. Bars represent the percentage of infected plants with 
symptoms indicated in the legend. N corresponds to the total 
number of plants infected. Deletion of the individual numl or kinl 
genes results in slightly reduced pathogenicity, whereas the numl I 



fa'«7-double deletion has a much stronger effect on pathogenicity. 
Depicted below are disease symptoms (third leaf below the site of 
injection) of representative maize plants seven days after 
inoculation. 
(PDF) 

Figure S19 Numl has neither influence on morphology or 
abundance of microtubules and actin elements, nor on the 
distribution of the Kinl motor protein. Microtubules were 
analyzed using a TubhmCherry fusion protein, the actin 
cytoskeleton was visualized using the lifeact-method (lifeactyfp) 
[44]. Both fusion proteins were expressed under control of the 
constitutively active P^promoter in strains AB31 and AB?AA- 
numl. For the localization of the Kinl motor protein in strains 
AB33 and AB'i'iAnuml, an N-terminal 3eGFP:Kinl fusion was 
used [14]. The fluorescent cytoskeletal elements were microscop- 
ically analyzed 12-14 hours after induction of the bEl/bW2- 
heterodimer in arabinose-containing CM-medium for AB31 and 
nitrate-containing minimal-medium for AB33, respectively. 
mCherry, YFP- and GFP-signals are shown. Scale bars: 10 |J,m. 
(PDF) 

Figure S20 Orientation of the microtubule cytoskeleton in AB33 
and AB33Anuml sporidia and hyphae. Orientation of the 
microtubule cytoskeleton was addressed in strain AB33PeblR_GT 
(AB33 pebl:mRFP potefGFP:Tubl) [14] and its Anuml derivative. In 
AB33PeblR_GT wild-type cells, the homologue of the microtu- 
bule (MT) plus-end marker EB1, Pebl, fused to RFP, localizes to 
growing MT plus-ends [45], which are predominantly oriented 
towards the cell poles of the hyphae [9,14], or, in interphase cells, 
towards the budding cell and the opposing cell pole of the mother 
cell [46]. (A) Motility of PebhRFP in budded cells of AB33- 
PeblR_GT (WT) and its Anuml derivative. PebhRFP binds to 
plus-ends of microtubules, and its motility (pictures were taken in 
15 sec intervals) indicates the orientation of the microtubule [46]. 
Arrowheads and asterisks indicate individual Pebl:RFP dots. Scale 
bar: 2 u.m (B) Quantification of microtubule orientation in 
interphase sporidia. AB33PeblR_GT (WT) and its Anuml 
derivative were grown in YEPSL medium to mid log phase. 
Black arrows and numbers above and below the cartoons give the 
percentage of Pebl -RFP dots moving to either tip of the cell. 
White arrows show the relative amount of Peb 1 -RFP to the bud 
neck. Numbers are based on 45 sec intervals of 30 and 46 cells in 
wild-type and /frami-strains, respectively. The total number of 
signals used for quantification and the number of cells is given in 
brackets. Quantification was done according to the protocol given 
in [46]. (C) In AB33PeblR_GT wild-type hyphae, Pebl-RFP is 
predominantly oriented towards the cell poles. (D) In AB33- 
PeblR_GT Anuml, Pebl:RFP labeled MT plus-ends are also 
found at septa (marked by arrowheads), which are often placed in 
the middle of the hyphae. In some hyphae, mitotic spindles were 
observed (marked by asterisks), implicating that the nwmi-deletion 
leads to a deregulated cell cycle arrest. Peb 1 -RFP and GFP-Tub 1 
fluorescence was analyzed 8-10 hours after induction of hyphal 
growth in nitrate-containing minimal medium. Scale bars: 10 u,m. 
(E) Schematic representation of the orientation of the microtubule 
cytoskeleton in wild-type and Anami-hyphae. Peb 1 -localization 
indicates microtubule plus-ends. 
(PDF) 

Figure S21 Influence of Numl on the localization of the 
myosin chitin synthase Mcsl. (A) In wild-type cells, 3eGFP:Mcsl 
localizes to growing bud tips, where it participates in the synthesis 
of chitin [15]. In contrast, A num 1 -mutants shows reduced Mcsl- 
accumulation at the growth region. Shown are DIC and GFP- 
images, as well as contrast-inverted images. Scale bars: 5 |J,m. 
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(B) Quantitative analysis of the 3eGFP:Mcsl signal intensity in 
growing bud cells. The bar diagram shows corrected total cell 
fluorescence (CTCF) of wild-type (control) and A numl -mutant 
cells. ImageJ software was used for quantification. N corresponds 
to the number of cells analyzed. (C) In AB33 hyphae, 
3eGFP:Mcsl localizes to distinct foci close to the cell membrane 
and forms a gradient towards the growth zone within the hyphal 
apex or localizes to the basal retraction septum (arrowhead), 
respectively. In contrast, in AB33Anuml, no tip-ward gradient is 
obvious and fewer foci are observed at the cell membrane. In 
many cases Mcsl does not localize to delocalized septa in the 
Anuml-mutant (lower left panel, arrowhead indicates a septum). 
Mcs 1 -fluorescence was analyzed 12 hours after induction of 
hyphal growth in nitrate-containing minimal medium. Scale bars: 
10 nm. (D) Enlarged images of the sections marked in (C). Images 
were contrast-inverted. 
(PDF) 

Figure S22 Numl influences the dynamic of the RNA-binding 
protein Rrm4. To analyze the impact of numl on Rrm4 
movement, an Rrm4:eGFP fusion protein was expressed in 
AB33 and AB33Anuml. bEl/bW2 were induced for 16 hours in 
nitrate minimal medium, and Rrm4:eGFP localization was 
analyzed by fluorescence microscopy. (A) In AB33, Rrm4:eGFP- 
fusion proteins localize on bidirectional moving particles that are 
evenly distributed within the hyphae, as previously described [46]; 
see Supporting Video S3. (B) In AB33Anuml the Rrm4-particles 
accumulate at the hyphal tip (asterisk, upper panel) or within the 
hyphae (lower panel); see also Supporting Video S4 Arrowhead 
indicates septum. Scale bars: 10 u,m. 
(PDF) 

Video SI Motility of Yupl-labeled early endosomes in AB33. 
Visualization of early endosomes in strain AB33 using a 
Yupl:eGFP-fusion protein under control of the constitutively 
active P^promoter as a marker. Depicted is the GFP-fluores- 
cence 14 h after induction of the bE/bW-heterodimer in nitrate- 
containing minimal medium. Rapid, bi-directional motility of 
early endosomes is visible throughout the filament. Immobile 
bright fluorescent spots most likely represent vacuoles, as Yupl 
was also reported to bind to vacuolar membranes [4]. The image 
sequence was taken using an Axioimager.Zl microscope equipped 
with an Axiocam MRm camera (Carl Zeiss). In total, 50 frames 
with an exposure time of 253 msec were taken every 0,46 sec for 
23,173 sec. Time is given in hrs:min:sec.msec. The scale bar 
represents 10 urn. 
(MOV) 

Video S2 Motility of Yupl-labeled early endosomes in AB33A- 
numl. Visualization of early endosomes in strain AW33Anuml using 
the YuplxGFP-fusion protein under control of the constitutively 
active P, )(? /-promoter as a marker. Depicted is the GFP-fluores- 
cence 14 h after induction of the bE/bW-heterodimer in nitrate- 
containing minimal medium. Endosomal trafficking is consider- 
ably reduced and vesicles accumulate at hyphal tips as well as 
around delocalized septae. The image sequence was taken using 
an Axioimager.Zl microscope equipped with an Axiocam MRm 
camera (Carl Zeiss). In total, 50 frames with an exposure time of 
44 msec were taken every 0,26 sec for 12,891 sec. Time is given in 
hrs:min:sec.msec. The scale bar represents 10 um. 
(MOV) 

Video S3 Motility of Rrm4-labeled particles in AB33. Visuali- 
zation of the RNA-binding protein Rrm4 in strain AB33 using an 
Rrm4:eGFP fusion protein that was expressed under control of its 
native promoter. Depicted is the GFP-fluorescence 16 h after 



induction of the bE/bW-heterodimer in nitrate-containing 
minimal medium. The protein localizes to rapid, bi-directionally 
transported particles. The image sequence was taken using an 
Axioimager.Zl microscope equipped with an Axiocam MRm 
camera (Carl Zeiss). In total, 20 frames with an exposure time of 
500 msec were taken every 0,96 sec for 19,266 sec. Time is given 
in hrs:min:sec.msec. The scale bar represents 10 urn. 
(MOV) 

Video S4 Motility of Rrm4-labeled particles in AB33Anuml . 
Visualization of the RNA-binding protein Rrm4 in strain 
AB33Anuml using an Rrm4:eGFP fusion protein that was 
expressed under control of its native promoter. Depicted is the 
GFP-fluorescence 16 h after induction of the bE/bW-heterodimer 
in nitrate-containing minimal medium. The particles are fre- 
quently observed in accumulations at various parts of the hyphae 
and trafficking is considerably reduced. The image sequence was 
taken using an Axioimager.Zl microscope equipped with an 
Axiocam MRm camera (Carl Zeiss). In total, 25 frames with an 
exposure time of 500 msec were taken every 1,16 sec for 
29,141 sec. Time is given in hrs:min:sec.msec. The scale bar 
represents 10 |a.m. 
(MOV) 

Text SI Prpl9:3xHA and Cefl:3xHA fusion proteins are 

functional in U. maydis. 

(DOCX) 

Table SI Intron retention in AB31 and AB3lAnuml. 
(XLSX) 

Table S2 Intron containing genes within the FunCats "cell 

cycle" or "DNA repair". 

(XLSX) 

Table S3 Splicing defect upon numl deletion in genes encoding 

transcription factors. 

(XLSX) 

Table S4 Expression of b-target genes in AB31 compared to 

AB3\Arbfl and AB3lAnuml. 

(XLSX) 

Table S5 Altered gene expression in AB31 Anuml . 
(XLSX) 

Table S6 Numl interacting proteins identified by Yeast Two- 
Hybrid screening. 
(DOCX) 

Table S7 U. maydis strains used in this study. 
(DOCX) 

Table S8 A. nidulans strains used in this study. 
(DOCX) 

Table S9 Oligonucleotides used in this study. 
(DOCX) 
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